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Abstract 

A systematic field-theory approach to Electron Spin Resonance (ESR) in the 
5 = 1/2 quantum antiferromagnetic chain at low temperature T (compared 
to the exchange coupling J) is developed. In particular, effects of a transverse 
staggered field h and an exchange anisotropy (including a dipolar interaction) 
5 on the ESR lineshape are discussed. In the lowest order of perturbation 
theory, the linewidth is given as oc Jh 2 /T 2 and oc (5/J) 2 T, respectively. In 
the case of a transverse staggered field, the perturbative expansion diverges 
at lower temperature; non-perturbative effects at very low temperature are 
discussed using exact results on the sine-Gordon field theory. We also compare 
our field-theory results with the predictions of Kubo-Tomita theory for the 
high-temperature regime, and discuss the crossover between the two regimes. 
It is argued that a naive application of the standard Kubo-Tomita theory to 
the Dzyaloshinskii-Moriya interaction gives an incorrect result. A rigorous 
and exact identity on the polarization dependence is derived for certain class 
of anisotropy, and compared with the field-theory results. 
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I. INTRODUCTION 



Quantum spin chains have been studied extensively for both their experimental and 
theoretical interests. Among many experimental methods of investigation, Electron Spin 
Resonance (ESR) is unique for its high sensitivity to anisotropy. While the theory of ESR 
has been studied 1-4 for a long time, there remain important open problems, especially for 
strongly interacting systems. One of the problems is that, generally one has to make a 
crucial assumption about the lineshape at some point during the calculation. As we will 
demonstrate, such an assumption could be incorrect in some cases although it might have 
been taken for granted in the literature. In addition, in an actual calculation one has 
to calculate various correlation functions. Traditionally, crude approximations such as the 
high-temperature approximation, the classical spin approximation and the decoupling of the 
correlation functions are used. However, these approximations break down when the many- 
body correlation effects are strong. As a consequence, rather little has been understood 
about ESR when many-body correlations become important. Even in the cases which were 
believed to be understood with the existing theories, there appear to be subtle problems. 

In this paper, we study ESR in S — 1/2 quantum spin chains in the "one-dimensional 
critical region" where the temperature T is sufficiently small compared to the characteristic 
energy of the exchange interaction J (but T is still large compared to three-dimensional 
ordering temperature or spin-Peierls transition temperature.) We stress that ESR in such 
a region is essentially a many-body problem. Here, many of the traditional theoretical 
techniques lose their validity. Instead, (1 + l)-dimensional field theory should describe the 
universal, low-energy/large-distance behavior. Our main purpose in the present paper is 
to develop a new approach to ESR based on field theory (bosonization) methods. At least 
for several simple cases (which are of experimental interest) we are able to formulate the 
problem in terms of the systematic Feynman-Dyson perturbation theory, avoiding previously 
made ad hoc assumptions. We will study several consequences of our theory for two types 
of perturbations of the one-dimensional S — 1/2 Heisenberg antiferromagnet: a staggered 
field, and an exchange anisotropy (or dipolar interaction). When the effect of the anisotropy 
is small, the ESR lineshape is shown to be Lorentzian up to a possible small smooth back- 
ground; the width and the shift of the Lorentzian peak are given perturbatively. In one 
dimension, it was argued that the diffusive spin dynamics leads to a non-Lorentzian line- 
shape, which is indeed observed in the S = 5/2 antiferromagnetic chain TMMC 4 . However, 
our results imply that the argument does not apply to the present case of the S — 1/2 
antiferromagnetic chain at low temperature. 

In a compound with a low crystal symmetry permitting a staggered component of the 
gyromagnetic tensor or a staggered Dzyaloshinskii-Moriya (DM) interaction, an effective 
staggered field is also produced by the applied uniform field. The staggered field corre- 
sponds to a relevant operator in the Renormalization Group sense, and is related to the 
field-induced gap phenomenon recently found in several quasi-one dimensional S — 1/2 
antiferromagnets 5 ' 6,8 ' 9 . Since it is a relevant operator, one may expect that its effect is en- 
hanced at lower temperatures. Indeed, we find that the staggered field contributes to the 
linewidth proportionally to h 2 /T 2 where h is the magnitude of the staggered field. We pro- 
pose this as an explanation of the peculiar low-temperature behavior 10 found in ESR on Cu 
Benzoate nearly 30 years ago. Moreover, we propose that the sharp resonance found at very 
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low temperature 11 , which was understood as a signature of a three-dimensional Neel order- 
ing, may well be understood in a purely one-dimensional framework based on sine-Gordon 
field theory. 

On the other hand, dipolar interactions or exchange anisotropies are present in virtually 
any real material. We find that their contribution to the linewidth is proportional to T, which 
appears to be consistent with existing experimental data on several quasi one dimensional 
S — 1/2 antiferromagnet such as CuGeO"3, KCUF3 and NaX^Os. 

Basic ideas and some of the results in the present paper were presented briefly in Ref. 12. 
This paper is organized as follows. In Section II, we briefly review the basics of ESR 
in interacting spin systems, including a few (apparently) new results, namely an exact and 
rigorous identity on the polarization dependence, and the relation between the Kubo-Tomita 1 
and the Mori-Kawasaki 2 ' 3 theories. In Sections III and V we develop a new framework for 
studying ESR in quantum spin chains, based on field theory methods and, in particular, the 
Dyson formula expressing the Green's function for a scalar field in terms of the self-energy. 
It is applied in Sections IV, VI, VII and VIII to systems with an exchange anisotropy (or 
dipolar interaction) or a transverse staggered field. (The case of an exchange anisotropy 
with the axis parallel to the field turns out to be easier to treat and not to require the self- 
energy formalism. Therefore it is treated first, in Section IV.) In Section IX, we compare 
our results to those in the high-temperature regime obtained with the previous approach. 
Section X is devoted to conclusions. Appendix A contains an alternative derivation of an 
old formula for the width/shift first derived by Mori and Kawasaki 2 ' 3 . 



II. ELECTRON SPIN RESONANCE 

A. Definition of the problem 

A single spin in a magnetic field H has energy levels separated by the Zeeman energy 
Ez = gusH- If an electromagnetic wave of angular frequency uj is applied to such a 
system, resonant absorption occurs when hu = Ez and the polarization (direction of the 
oscillating magnetic field) is perpendicular to the static field. When the spins are coupled 
by interactions, the physics is of course not that simple. However, generally some resonant 
absorption occurs also in the interacting system. This is the phenomenon of ESR which we 
study in the present paper. In an interacting system, it is also possible to observe absorption 
of the electromagnetic wave polarized parallel to the static magnetic field (so called Voigt 
configuration.) In this paper, we focus on the standard (Faraday) configuration, which 
measures the absorption of the electromagnetic wave polarized perpendicular to the static 
magnetic field. 

Assuming that the absorption can be described by linear response theory, the absorption 
intensity I{uj) per volume for the radiation linearly polarized in the a _L z direction is given 
by 

/H = ^xL(g = o,^), (2.1) 

where Hr is the amplitude of the radiation and x" is the imaginary part of the dynamical 
magnetic susceptibility, x" is related to the retarded Green's function Q R as 
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(2.2) 



where is defined by 

g%,(q,u) = -i / dtJ2([S a (x,t),S^0,0)])e- iqM (2.3) 

where (. . .) is the statistical average at temperature T. In most experiments, the applied 
electromagnetic wave is typically in the microwave regime, and its wavelength is very large 
compared to all relevant length scales in the antiferromagnet since the spin-wave velocity is 
much less than the speed of light. Thus, in ESR, the dynamical susceptibility is measured 
at essentially zero momentum q = 0. ESR probes the dynamics of the system only at 
the special momentum q = 0, in contrast to neutron scattering which can be used to scan 
momentum space. However, as we will explain below, there is an interesting feature at the 
special momentum q = 0. Together with the relatively easy availability of highly precise 
data, ESR offers a unique insight into magnetic systems which would be difficult to obtain 
with other experimental methods. 

A remarkable feature of ESR is that, if the Hamiltonian of the system (apart from the 
Zeeman term) is isotropic (ie. £77(2) symmetric), the resonance is still at the Zeeman energy 
and completely sharp, as if there is no interaction at all. This result can be deduced rather 
easily from the equation of motion, as we will show in the following. Throughout this paper, 
we take the direction of the static applied field as the z-axis. Let us consider the total 
Hamiltonian 



h = h + h z , (2.4) 

where Hz = —HJ^jSj is the Zeeman term and Ho is the exchange Hamiltonian which 
is assumed to be £77(2) symmetric. We choose units so that h = gfis — 1 except where 
explicitly mentioned otherwise; these constants can be recovered by dimensional analysis. 
We also assume H is much smaller than the exchange coupling J (ie. energy scale of 
Ho). As we have mentioned above, in ESR the electromagnetic wave is coupled to the 
q = component of the spin operators, namely the total spin operators £ Q = J2j The 
Heisenberg equation of motion for S + = S x + iS v reads 

7S*+ 

__ = i[n, £+] = i[H z , S + ] = -iHS + , (2.5) 

because Ho commutes with £ + due to the £77(2) symmetry of Ho. It follows that 
S + (t) = S + e~ tHt , and consequently x + ~(0,cj) oc 5(uj — H). This means that the resonance is 
completely sharp, and located exactly at the Zeeman energy. Namely, this resonance has the 
lineshape identical to ESR in a single (non-interacting) spin in spite of an arbitrary strong 
exchange interaction. On the other hand, the absorption intensity is generally affected by 
the exchange interaction Ho- For example, in a spin-gap system at zero temperature, the 
absorption intensity is zero if the applied field H is smaller than the gap. 

As we have seen, the completely sharp resonance is related to the £77(2) symmetry of the 
exchange Hamiltonian Ho- While it is natural that symmetries of the system are important 
in determining the dynamics of the system, the present situation is rather unique, for the 
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577(2) symmetry is explicitly broken down to U(l) by the applied static field but is still 
essential in ESR. This peculiar feature is related to the fact that the applied field couples to 
the total magnetization S z = J2j 5|, which is a generator of the global SU (2) symmetry and 
is conserved under Tio. Since the total magnetization and Hamiltonian are simultaneously 
diagonalizable, the applied field does not change the eigenstates of the system, if they are 
classified by S z . The only effect of the static applied field is to shift the energy levels of the 
eigenstates; the shifted energy levels still reflects the 577(2) multiplet structure. This kind 
of "weak" symmetry breaking by one of the symmetry generators preserves some structures 
of the fully symmetric system. In ESR of an isotropic system Tio, the 577(2) symmetry is 
only weakly broken and is essential in determining the ESR spectrum. 

A similar application of the concept of weakly broken global symmetry was also exploited 
recently by S. C. Zhang 13 in his 50(5) theory of high-T c superconductivity. Namely, in the 
50(5) theory, the most important terms in the effective Hamiltonian are 50(5) symmetric 
one and the chemical potential couples to one of the generators of the global 50(5) symmetry. 
The so-called 7r-excitation in this context is a sharp resonance which is similar to ESR in 
isotropic spin systems. 

In real magnetic systems, there are various types of anisotropy, such as the dipolar 
interaction. Let us write the total Hamiltonian as 

n = n + n' + n z , (2.6) 

where Ti' is the symmetry-breaking perturbation. Throughout this paper, we assume the 
interaction to be nearly isotropic, namely that Ti' is small compared to the other terms Tio 
and Tiz- Once the perturbation Ti' is added, the argument leading to the delta-function 
resonance at the Zeeman energy breaks down. Thus, in general, the addition of Ti' causes 
changes in the lineshape, such as a broadening and a shift of the resonance. The main 
theoretical problem is then to calculate the absorption spectrum for the given Hamiltonian 
Ti and other conditions such as the temperature of the system. 

B. Previous theories 

The existing approaches to ESR, such as those of Kubo and Tomita 1 and of Mori and 
Kawasaki 2 ' 3 were developed mainly during 1950s-60s. Here we summarize briefly, a part of 
those achievements which is closely related to our analysis. 

When the isotropic exchange interactions between spins are weak, namely Tio is much 
smaller than Tiz, the lineshape is generally expected to be Gaussian. On the other hand, 
once the anisotropy Ti' is present, strong isotropic exchange interactions Tio between spins 
affects the ESR lineshape, even though it does not break the SU (2) symmetry by itself. In 
the presence of the strong interaction (Tio ^ Hz), which applies to the problem considered 
in this paper, the lineshape is generally expected to be Lorentzian. (On the other hand, 
the lineshape is believed to be neither Gaussian nor Lorentzian, when the spin diffusion 
is dominant. 4 ) The effect of the isotropic exchange interactions on the lineshape has been 
traditionally called Exchange Narrowing. We emphasize that ESR in such an interacting 
spin system probes the collective motion of the many-body system. In this paper, we focus 
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on this limit of strong isotropic exchange interaction, while other cases have been discussed 
previously as well 1-3 . 

For the case of the Lorentzian lineshape, Mori and Kawasaki 2 proposed a formula, which 
we call the MK formula, for the linewidth rj: 

V = T^-ErM-GHuiu = H)], (2.7) 

where Xu is the magnetic susceptibility and G^ Ai (u) is the Fourier transform of the unper- 
turbed retarded Green's function 

G^(t) = - l 9(t)([A(t),AHO)])o, (2.8) 

where (. . .) is the expectation value under the unperturbed Hamiltonian 7i + 'Hz, 0(t) is 
the step function, and A is defined by the commutator 

A=[H',S + ]. (2.9) 

In this paper, Q refers to a full Green's function calculated using the Hamiltonian including 
the perturbation Ti', while G denotes the unperturbed Green's function evaluated in the 
absence of the perturbation. Both kinds of Green's functions (Q and G) should be evaluated 
including the Zeeman term Hz, in the original spin chain context. However, as we will see in 
Section III, in the effective field theory, the Zeeman term is absorbed by a momentum shift. 
Thus, the Green's functions in the effective field theory will be defined without explicitly 
including the Zeeman term. 

In addition to the broadening, the perturbation Ti! also causes a shift of the resonance 
energy; the shift is given by 

Au = ^=L { {[A , S-]) ~ ReG^(u = H)}. (2.10) 

This formula for the shift is slightly different from the one given in the original paper 2 . We 
believe that ours is the correct one in the lowest order of perturbation theory. 

The derivation of the MK formulae in the original paper seems somewhat involved, and 
it is not clear to us what assumptions are necessary to prove them. However, we found 
that the MK formulae is indeed exact in the lowest order of the perturbation theory, if the 
(single) Lorentzian lineshape is assumed. Explicitly speaking, we must assume: 

where £ is a smooth function of uo near u = H. Regarding £ as a constant near the 
resonance, we obtain a Lorentzian line-shape. [Setting £ = in Eq. (2.11) gives the exact 
result for the isotropic case TC = 0.] The simple, and possibly new, alternative derivation 
using the equation of motion is presented in Appendix A. 

On the other hand, Kubo and Tomita 1 studied ESR using a somewhat different formu- 
lation. For the case of Lorentzian lineshape, their theory gives the following formula for the 
linewidth, at high temperature: 
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where the expectation value is the static correlation function. We shall call this the KT 
formula in this paper. We could not find in the literature how the two formulae (2.7) and 
(2.12) are related. On the other hand, if the KT formula (2.12) for the Lorentzian lineshape 
is indeed valid at high temperature, it must be consistent with the MK formula. In fact, 
we have verified that the KT formula (2.12) follows from the high-temperature limit of the 
MK formula (2.7) with a certain assumption. The derivation is given as follows. Taking the 
Fourier transform of eq. (2.8), at temperature T, 

GltfM = ~ f dte^Tr([A(t),AHO)]e-^ +n ^ T ), (2.13) 
where Z = Tre~^ ()+Wz ^ T . Expanding this up to the first order in 1/T, we find, 
GiitM ~ ^- f dm([H + n z ,A(t)]A\0))e^ 

dm{^{t)A\^e wt 

TZ Jr [A(0)A\0)} - i£- jH dtTr [A(t)Am] e**, (2-14) 

where the time evolution is defined with respect to the unperturbed Hamiltonian 7i + 'Hz 
and — Trl is the partition function in the infinite temperature limit. The first term 
is real and does not contribute to the imaginary part. If we assume that the dynamical 
correlation function at infinite temperature Tr ^(t)^4(0)J decays exponentially with the 

characteristic time r c ~ 1/J, the second term gives — i-jf(A^(0)A(0))oo, where ()oo is the 
expectation value at the infinite temperature and we use u <C J . We note that a similar 
assumption was made also in the original derivation of eq. (2.12) in the Kubo-Tomita paper. 
Thus the MK formula (2.7) reduces, in the high-temperature limit, to 

Because (S + S~) ~ 2% U T in the high-temperature limit, this is equivalent to the KT for- 
mula (2.12). We note that, because r c ~ 1/ J is valid only as an order-of-magnitude estimate 
at best, the KT formula has the uncertainty of an overall constant factor. 

Recently, a numerical approach to ESR in quantum spin chains is also being developed 14 
by a direct calculation of the dynamical susceptibility Since it is based on an exact 

diagonalization of the full spectrum of short chains, it is restricted to rather short chain of up 
to 10 spins even for S — 1/2, making finite size effects rather severe. On the other hand, the 
direct numerical calculation is applicable at any temperature. In contrast, the field theory 
approach, which we will develop in the present paper, is valid only at low temperatures while 
it is based on the thermodynamic limit. Thus they are complementary to each other. 

We remark that some results quite closely related to ours were derived by Giamarchi and 
Millis 30,31 in their work on the ac conductivity of a Tomonaga-Luttinger (TL) liquid. We 
will comment on the connections with our work later. 
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C. Polarization dependence 



When observing ESR in the Faraday configuration, the polarization of the electromag- 
netic wave is perpendicular to the direction of the static magnetic field, which we take as 
the z-axis. There are still two independent possible polarizations; the linear polarization 
can take any direction in the xy-plane. Except when the total Hamiltonian 7i is invari- 
ant under a rotation about the z axis, the absorption spectrum generally depends on the 
polarization. Within the linear response theory, the dependence comes from the difference 
between the dynamical susceptibility x"x(0,w) ^ x"y y (Q,uo). The MK formula ignores the 
possible polarization dependence, because it deals with x%- ~ Xxx + Xyy, an d not x'xx an d 
Xyy separately. The polarization dependence was discussed theoretically first by Natsume et 
al. 15-17 generalizing the Kubo-Tomita theory. It has also been observed experimentally 15,16 
and numerically 14 . 

However, apparently it has been not recognized that, for some special cases, an exact 
and rigorous result on the polarization dependence can be derived easily from the equation 
of motion. Let us consider the special case in which the perturbation 7i' is written in terms 
of the x component of the spin operator SJ. The examples include the transverse staggered 
field H! = h^2j(—iySj, and the exchange anisotropy with the anisotropy axis in the x 
direction W = 5J2j Sj&j+v 01 these cases, [S X ,H'} = holds, and consequently 



This identity leads to 



d^ x 

= HS y . (2.16) 

dt 



xL(0»w) = ^xi(0,o;), (2.17) 



and more generally, for the polarization in the direction a in the xy-plane, 

rt^^w^rt,,, (2 , 8) 

where $ is the angle between x and a. (In the notation of Refs. 15-17, 6 = 90 degree and 
their corresponds to our $.) 

For a sharp resonance concentrated near uj ~ H, the polarization dependence is not 
significant. However, if the center of the resonance is defined by the average frequency of 
the spectrum 

f ujl a (uj)duj 

a« = T7TTi ' 2 - 19 

J I a (uj)duj 

oj x < oj y because the higher frequency part is emphasized in the latter. As a consequence, 
there is a positive frequency shift for the polarization in y axis, compared to the case where 
in x axis. 

This is in agreement with theoretical and experimental results in Refs. 15, 16 and numer- 
ical results in Ref. 14, on the exchange anisotropy. We note that, in the actual experiments 
on ESR, the resonance frequency is kept fixed and the applied field is scanned to measure 
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the absorption. Because of this, it is customary to discuss the shift of the resonance field for 
a fixed frequency. The direction (positive or negative) of the field shift is opposite to that 
of the frequency shift we discuss in this paper. They find that the resonance field is shifted 
negatively for the polarization in y direction compared to the x polarization case, which 
is indeed consistent with our result. The angular dependence is also consistent with the 
theoretical formula in Ref. 15. On the other hand, in Ref. 17, the polarization dependence 
is studied by a different formalism (Mori's memory function method.) When the anisotropy 
axis is perpendicular to the applied field, the obtained polarization dependence is rather 
opposite to the above, and is in contradiction to our rigorous result (2.17). 



III. FIELD-THEORY APPROACH TO THE 5=1/2 HEISENBERG 
ANTIFERROMAGNETIC CHAIN 



A. Bosonization of S = 1/2 Heisenberg chain 



In the present paper, we mainly discuss ESR on the one- dimensional S — 1/2 Heisenberg 
ant if erromagnet 

Ho = Sj ■ S j+ i. (3.1) 

3 

with symmetry-breaking perturbation H! and of course the Zeeman term Tiz- The low- 
energy physics of the one-dimensional quantum antiferromagnets is well described by field 
theory methods (bosonization). In this section, we briefly summarize the aspects of this 
approach that are relevant to the present discussion of the ESR. We refer the reader to 
Refs. 18, 19 for more details. While the method is now standard, here we also clarify 
subtleties specific to ESR problems, which are related to the weakly broken SU (2) symmetry 
discussed in Section II. 

The effective field theory of the S — 1/2 Heisenberg chain 7i is given by the free boson 
Lagrangian 

C = \ [(Mf - (<9i0) 2 ] , (3.2) 

where x° = vt, x l = x and we make identification <fi ~ <fi + 2irR with the compactification 
radius R. The radius R is actually fixed to the value l/\^2n by the SU(2) symmetry. 
Hereafter we set v — 1 for simplicity; the spinon velocity v can be recovered by dimensional 
analysis when necessary. 

At zero uniform field, the spin operators may be written in terms of the field as follows: 

SJ ~ C-e- i2wR * cos | + C;(-iy e - t27TR *, (3.4) 

where the dual field is defined in terms of right-mover and tp^ as = ipn + LpL and 
— ¥r ~~ Vl- While S z and S x ' y are represented in a very different way, their correlation 
functions turn out to be equal at the SU(2) invariant radius R = l/y/2n, as required from 
the symmetry of the original Heisenberg chain. 
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FIGURES 




FIG. 1. The spin structure factor of the S = 1/2 Heisenberg antiferromagnetic chain at T = 0. 
It is non- vanishing only in the filled region shown in the frequency-momentum plane. The structure 
factor becomes a delta function 8{u — vq) in the q — > limit. 



The dynamical structure factor S aa (Fourier transformation of the spin correlation func- 
tion) of the Heisenberg chain has been studied in detail. It is equivalent to the dynamical 
susceptibility for T = and u > 0. At zero temperature, the dynamical structure factor is 
non-vanishing only in the limited region of the frequency uj - momentum q space shown in 
Fig. 1. The field theory actually can handle only the low-energy excitations near momen- 
tum and 7r. The structure factor for S zz near q = and q = n is given by the correlation 
function of d(p/dx and cos (<f>/R) respectively. At T = 0, they read 



S zz (uj,q) oc 6(u- \q\)6(u) (3.5) 



for q ~ and 



S*>, ?) oc . 6{u -\q-n\) (3.6) 

for q ~ 7r. It is noted that the structure factor is completely sharp and is delta-function like 
at q ~ 0. In fact, the structure factor at q ~ remains so even at finite temperature. As 
mentioned before, the structure factor is of course isotropic (S xx = S yv = S zz ) at H = for 
the isotropic Heisenberg chain. 

Now let us consider the effect of the applied magnetic field. The Zeeman term Hz in the 
Lagrangian becomes, upon bosonization, 



This term can be eliminated by a redefinition of the boson field 



4>{t,x)^4>{t,x) + -j=x. (3.8) 



10 



but remains unchanged. This is equivalent to the shift of chiral fields as 

<PR^<PR + -z^j=Hx,ip L ^<p L + —=Hx. (3.9) 



While this leaves the free Lagrangian unchanged, it does change the bosonization formulae 
of physical spin operators: 

ss ~ m+ 2^S + c i cos| l + (i/ + ^- (310) 

S ± ~ C'e-'^cos + Hx) + C;(-l) j e~^ R t (3.11) 

The first term m in S z represents the expectation value of the magnetization induced by 
the magnetic field H. For a small magnetic field, m is proportional to the field H . Another 
important feature is that the applied field induces the shift of the soft-mode momentum 20 ' 21 . 
The shift occur differently for the longitudinal (z) and the transverse (x, y) components. 
The gapless points under the applied uniform field H are at q = (uniform part) and 
q = 7r ± H ( "staggered" part) for the longitudinal modes. For the transverse modes, they 
are at q — ±77 ("uniform" part) and q = -n (staggered part.) 

Let us focus on the transverse mode near q — 0, because the transverse mode at q = is 
measured in ESR in the Faraday configuration. For simplicity, here we restrict ourselves to 
zero temperature. In the low energy effective theory, the "uniform" part of the 5 ,± is given 



where we have used the £77(2) symmetric compactification radius R = 1/V27T (see below for 
reason for taking this value.) This gives the correlation function of S ± at zero temperature: 

(S+(t ,,) 5 -(0,0))« I? -^ + I ^!_ (3 . 13 ) 
Dynamical structure factor £+_, which is the Fourier transform of the above is, 

AHx „—iHx 



/oo roo 
dt / dxe l{ult - qx) 
-oo J —oo 

/oo 
dx[-2nu9(cu)]e" tqx 
-oo 



e 

+ 



(t — ie + x) 2 (t — ie — x) 5 

iHx~iuix _|_ ^—iHx+iujx 



oc u6(uj)[5(uj - H + q) + 5(lo - H - q)]. (3.14) 

The other one £_+ is given by replacing H — > —H in the above, using the time reversal 
transformation. Thus 

S-+(u,q)oLu6(u)[8(u + H + q) + 8{u) + H-q)\. (3.15) 

Namely, £ + _ and £ h give different branches of excitation. The fact that £_ + does not 

contain the branch (3.14) was recognized earlier (see Fig. 17 of Ref. 21.) On the other 
hand, that £ + _ lacks the branch (3.15) (at least in the low-energy limit) was apparently not 
appreciated in Fig. 18 of Ref. 21. S xx and S yy are given by their superposition 
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S xx (u, q) = S yy (uj, q) oa u[S(u - \q + H\) + S(u - \q - H\)]. 



(3.16) 



This zero-temperature transverse structure factor near q = under the applied magnetic 
field is shown in Fig. 2. Because the structure factor near q = was sharp, and the gapless 
point is shifted by H, we expect a sharp resonance at energy uo ~ H at q = 0. This 
corresponds to the expected paramagnetic ESR for the isotropic Heisenberg chain. 



FIG. 2. The zero temperature transverse spin structure factor S xx (uj,q) = S yy (uj,q) of the 
S = 1/2 Heisenberg antiferromagnetic chain under an applied field H, near q = 0. It is approxi- 
mately proportional to lo[S(uj — \q — H\) + 5{u — \q + H\)\, giving the resonance at q = 0, uo = H. 
This consists of two branches coming from S + - and S-+, which are marked by H — and — h in the 
graph. In fact, there is a small spreading of the spectrum and the structure factor is generally not 
a perfect delta function. However, it is exactly the delta function 5(u> — H) at q = 0, as explained 
in the text. 



However, it should be noted that we have so far ignored various renormalization effects 
due to the applied magnetic field. There are irrelevant operators, which themselves vanish 
in the low-energy limit but renormalize parameters of the low-energy effective theory. The 
way they renormalize is affected by the applied magnetic field. In general, the precise value 
of the momentum shift is given by 2nm rather than H, where m is the magnetization. This 
can be derived from the shift of Fermi momentum in the Jordan- Wigner transformation, and 
also is required from a rigorous version of Luttinger's theorem in one dimension 22 . Restoring 
the spinon velocity v, the ESR frequency appears to be given by 2rrmv. For the standard 
Heisenberg antiferromagnetic chain in an applied field, the magnetization m and the spinon 
velocity v can be obtained as a function of H from the Bethe Ansatz integral equation. 
Generally, 2nmv is different from H except in the zero field limit, implying that the ESR 
frequency deviates from H. However, this cannot be true, because the equation of motion for 
the original Heisenberg model (under an applied field) requires the resonance to be exactly 
at the frequency H . The resolution is that, the dispersion relation for q ~ is not completely 
linear. The curvature of the dispersion comes form irrelevant operators which break Lorentz 
invariance. Because of the curvature, the resonance frequency at q = is modified from 
2nmv, which is derived assuming the linear dispersion. What the equation of motion tells 
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us is that these renormalization effects miraculously cancel, to give the resonance exactly 
at uj = H for q = 0. With this nontrivial mechanism in mind, we will take the momentum 
shift as H, setting the spinon velocity v — 1. 

There is another "miraculous" cancellation similar to the above. At zero field, the 
compactification radius of the effective field theory is fixed to the special SU (2) symmetric 
value R = l/y/2n, as is required from the ££7(2) symmetry of the original Heisenberg 
model. However, in the presence of the applied field, the £77(2) symmetry is of course 
broken down to £7(1). Correspondingly, the radius R is renormalized away from the ££7(2) 
point by the applied field. The renormalized radius R as a function of the applied field H 
has been also obtained from the exact Bethe Ansatz solution 23 . It is indeed rather sensitive 
to H for small H/ J. A consequence of the radius renormalization is the dependence of the 
correlation exponents on the applied field. In particular, the "uniform" part of the transverse 
spin operator, which is relevant for ESR, is represented by the vertex operator of the type 
exp [±2niR(p ± <j)/R\\ its conformal weight is given by 



(A, A) = (1 + A', A') (3.17) 



or (A', 1 + A'), where 



A . = (frfl-'/fl)', (3 . 18 ) 

8n 

which does depend on R. As a result, the structure factor is no longer given by a delta- 
function for R ^ l/v^2~7r. More explicitly, the retarded Green's function of a conformal 
primary field with conformal weight (A, A) at finite temperature T is obtained explicitly 24 

as 

Gf A> A)(", Q) = ~ (2 7 rA)(2 7 rT) 2 ( A+A - 1 )5(A - i^±?, 1 - 2A)5(A - i^-, 1 - 2A), 

(3.19) 

where B denotes the Euler Beta function: 

and T is Euler's Gamma function. Considering the momentum shift induced by the applied 
field, the absorption measured in ESR corresponds to the Green's function evaluated at 
q = H. Thus, the spectrum is given by the delta function only if (A, A) = (1,0) or (0, 1), 
namely R = l/y/2n. The renormalization of R due to the applied field seems to imply 
that the ESR spectrum should not be given by a delta-function, even in the absence of the 
perturbation Ti! . 

However, this is inconsistent with the equation of motion of the original Heisenberg 
model. It predicts a completely sharp (delta-function) resonance precisely at the Zeeman 
energy even for a finite field H . Since the equation of motion is exact and rigorous for 
the original spin problem, we conclude that we should take the unrenormalized, ££7(2) 
symmetric value R = even in a finite field, for the calculation of the ESR. This 
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appears contradictory to the well-established renormalization of R due to the applied field. 
This is not a real contradiction, however, because the standard result on the renormalization 
of the radius is determined at the zero energy limit, while the ESR probes the excitation 
at the finite energy H. In general, effective coupling constants depend on the energy scale 
as a consequence of the renormalization. We may introduce an effective radius R{uS) as a 
function of the energy scale u. While the determination of the function R(u) in general is 
a tedious task, the exact equation of motion on ESR gives the restriction at the Zeeman 
energy: R{u = H) = 1/ \/2tx. The non-renormalization could be related to the qualitative 
understanding of the RG flow in the presence of the applied field, Fig. 7 in Ref. 7. The RG 
flow in the presence of the applied field is almost identical to that in the zero field, down to 
energy scale of 0(H), where the flow is "cut off." If we look at the energy H, the effective 
theory may be almost identical to the isotropic one. This argument would not, however, 
explain why the effective radius should be exactly at the SU(2) point. From the viewpoint 
of the field theory, this is again a miraculous cancellation between the renormalization by the 
uniform field and that by the finite energy. The equation of motion, although quite simple, 
gives an exact and highly nontrivial constraint on the effective field theory description. 

Thus, in the following calculations we do not include the radius renormalization due 
to the applied field, and take the 5'?7(2)-symmetric value R = l/\f2n. As a result, the 
appropriate effective field theory of ESR is an SU(2) symmetric one, namely the level- 
1 SU(2) Wess-Zumino-Witten (WZW) theory, even in a finite field; all the effects of the 
applied field are represented by the shift of the <ft field (3.8), resulting in the momentum 
shift (3.10) and (3.11). This may be regarded as a field theory representation of the crucial 
SU (2) symmetry which is broken only weakly, discussed in Section II. 

It is often convenient to introduce the operators in non-Abelian bosonization to make 
the symmetry manifest. SU (2) current operators J a (a = x, y, z) are related to the Abelian 
bosonization as follows: 

J' R {w) = iV^dip R {w), (3.21) 
J±(w) = ^gWsW^ ( 3>22 ) 
J*(w) = -iV&rdip L (w), (3.23) 
J$(w) = ^/2e TlV ^^\ (3.24) 

where J ± = J x ± iJ y , we have introduced complex coordinates w = r + ix (r = it) and 
(p(w,w) = (f(w) + (p(w) . Jr^l) is the right-mover (left-mover) component of the current, and 
we have normalized them by 

(r R ( Wl )J b R (w 2 )) = {wi _ W2)2 , (3-25) 

where a,b = x,y, z and the complex coordinate w = r + ix = —i(t — x) and likewise for the 
L sector. (We note that this is different normalization from Ref. 25.) 

The "uniform" part of the spin operators S a correspond to the SU(2) currents J a , while 
the "staggered" part is related to the SU (2) triplet n a = Trga a where the 577(2) matrix field 
gp(x, t) is the fundamental field of the Wess-Zumino-Witten non-linear cr-model. Eqs. (3.10) 
and (3.11) may be rewritten as 
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S z ~ -=( J Z R + J Z L ) + C s [cos (if + n)xn z + sin (ff + Tt)xtrg\ , (3.26) 

V8"7T 2 

S± ~ "7=^(4^ + Jl^)- + ("im^, (3.27) 

The "staggered" part of £ z may be written as (— l) x n z at H = 0, but is a mixture of n z and 
trg in a finite field. 

The ESR absorption intensity is related to the Green's function of S x,y ; thus what is 
needed in the field theory is the Green's function of J x,v at momentum iff. 



B. Perturbations 

Having established the effective field theory for the unperturbed system 7f o + 'Hz, we 
now want to calculate the effects of the perturbation H! on the ESR lineshape. Assuming 
that the perturbation Ti' is small, H! can be mapped to an operator of the level-1 SU{2) 
WZW theory. 

In principle, an infinite variety of symmetry breaking perturbations Ti' is possible. In 
fact, there are infinitely many operators also in the field theory. However, most of the 
operators have large scaling dimensions, and thus renormalize rapidly to zero under the RG 
transformation. Thus, at low enough temperatures, only a few types of perturbations with 
smaller scaling dimensions are important. 

The operators with the lowest scaling dimension 1/2 are n a and tig in WZW theory. 
In the original spin chain Hamiltonian (at ff = 0), they correspond to the staggered field 
(3 independent perturbations corresponding to three directions) and the bond alternation. 
However, the bond-alternation does not break the SU(2) symmetry and hence should not 
affect the ESR lineshape, although it is not trivial to see this in the field theory. On the 
other hand, the staggered field perturbation does break the £77(2) symmetry and thus affects 
the ESR lineshape. The operators of interest with the second lowest scaling dimension 2, 
which are marginal, are Jf;J R - They correspond to the exchange anisotropy in the spin chain 
Hamiltonian. We will discuss these two most important cases in later sections. 

While we use the SU (2) symmetric field theory, care should be taken with the momentum 
shift due to the applied field. The momentum shift is determined by a simple rule in Abelian 
bosonization formulation (3.9). Namely, if one writes some operator at zero field in terms 
of (/7s, the above replacement gives a correct formula under the finite field ff . The operator 
corresponding to the perturbation Ti! may contain an oscillating factor. While such a term 
may be ignored in order to know whether there is a finite excitation gap above the ground 
state, it should be retained in theory of ESR which probes finite momentum of the effective 
field theory. For a general perturbation, the oscillating factor appears in the effective field 
theory, and it makes the theoretical analysis rather complicated. In this paper, we focus on a 
few simple cases in which there is no oscillating term (with finite momentum) in the effective 
Lagrangian. This still includes several cases of physical interest which are mentioned below. 
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1. Transverse staggered field 



A quasi one-dimensional spin system often has an alternating crystal structure along 
the chain. In such a case, generally we expect two features which are absent in a uniform 
system. 

staggered g-tensor The magnetic field H couples to the spin as ^BJ2j, a ,b^ a [9ab + 
(—iy g s ab ]S b p where g s is the staggered component of the (/-tensor. 

Dzyaloshinskii-Moriya (DM) interaction The low symmetry allows the antisymmetric 
interaction 26 ' 27 J2j Dj • (Sj x S j+1 ) . 

The DM interaction can be either uniform (Dj = D) or staggered (Dj = (—1) 3 D.) 

When the staggered (/-tensor is present, an effective staggered field oc g s H is produced 
upon an application of the external field. The direction of the staggered field is often 
approximately perpendicular to the applied field, although it is not necessarily so. The 
effect of the DM interaction is less trivial, but it can be actually eliminated by an exact 
transformation. Let us consider the case of a staggered DM interaction, and choose the axes 
so that the DM vector D is parallel to the z-axis. Then the Hamiltonian including the DM 
interaction is given by 

n DM = jJ2Sj ■ s j+1 + (-iyD(s*s] +1 - s]s* +1 ) 

3 

= \ Et^-i^ + rstjS- w + (h.c.)] 

3 

+ J Yl\- S 2j-l S 2j + S 2j S 2j+l}^ (3.28) 
3 

where J = J + iD. Now let us define the angle a = tan -1 Dj J, and rotate the spin at site 
j by the angle (— l) J, o;/2 about z axis: 

S+ -> S+e i{ -^ a / 2 (3.29) 

Then we obtain the Hamiltonian of the XXZ chain 

H = Y\JS-S Z j + i + l -Q(S+ST +1 + h.c.)}. (3.30) 

3 ^ 

It is argued 28 that this anisotropic exchange can cancel the pre-existing one. 

Now suppose that an external field H is applied in x direction. The applied field is 
transformed as 

- £[cos| S* + (-iy sm^S]}. (3.31) 

j j 

by the above transformation. Thus, in the presence of the Dzyaloshinskii-Moriya interaction, 
the applied uniform field produces an effective staggered field 6 . For general orientations of 

— * 

D of the staggered DM interaction 



16 



H DM = '£(-iyD-(S j xS j+1 ), 

j 



(3.32) 



the effective staggered field due to the DM interaction is given by D x H/(2J) 

These two effects give an effective transverse staggered field which is approximately 
perpendicular to the applied field. This mechanism is important in studying properties 
of several quasi-one dimensional antiferromagnets including Cu benzoate 5-7 , Yb4As3 8 and 
Pyrimidine Cu dinitrate 9 . 

2. Exchange Anisotropy 

The exchange anisotropy is the second relevant perturbation which affects the ESR line- 
shape. The dipolar interaction which exists in any real magnetic system is given by, restoring 
the Bohr magneton /ib, 

H dp = {g^ B ) 2 Yl 

ij 

where represents the vector from site i and j and for the simplicity the ^-factor is assumed 
to be uniform and isotropic. In a spin chain, the vector is parallel to the chain direction, 
and the dipolar interaction reduces to an effective exchange anisotropy parallel to the chain 
direction. The effect would be essentially the same with the nearest-neighbor anisotropic 
exchange interaction, because the dipolar interaction strength decreases rapidly with the 
distance. 

Let us consider the simplest case of the exchange anisotropy 

H a = 5j2Sl l SJ +1 (3.34) 

j 

with a symmetry axis n, which effectively covers the case of the dipolar interaction if n 
is taken to be the chain direction. Even in this simple case, a variety of configurations is 
possible by changing the relative direction of n and the direction z of the applied field, as is 
often done in experiments. 

As mentioned before, for a general direction, the perturbation in the field theory is rather 
complicated, making a calculation from first principles difficult. Thus, in this paper, we will 
focus on the two simplest cases, namely when n || z and n -L z. The case n \\ z allows us a 
direct calculation of the lineshape and will be discussed in Section IV. The latter case n _L z 
will be discussed in Section VI, based on the self-energy approach developed in Section V. 

IV. EXCHANGE ANISOTROPY PARALLEL TO THE FIELD: DIRECT 

CALCULATION 

Here we consider the case where the anisotropy axis is parallel to the applied magnetic 
field, namely n = z in (3.34). In this case, it is obvious that there is no polarization 
dependence as S x and S y are equivalent. 



— * — * — * — * 
Si ■ Sj 3(S*i • ^ij)(>Sj 



(3.33) 
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In this case, the perturbation in the effective field theory is given, at zero magnetic field, 

as 



£a = -\J Z R Jl (4.1) 

where A is a parameter proportional to 5/ J, for a small anisotropy 5/ J. The proportionality 
constant A J/ 5 is non-universal and model-dependent. (For the standard Heisenberg antifer- 
romagnetic chain, A is determined in Section VIC together with a logarithmic correction.) 

Before performing an explicit calculation, let us see what can be said about the temper- 
ature dependence of the linewidth from a general scaling argument. The perturbation (4.1) 
is a marginal one with the scaling dimension 2. Thus, ignoring the logarithmic corrections, 
scaling arguments imply that the linewidth takes the scaling form 

V = Tf(j,j), (4.2) 

where we have used the fact that 77 has the dimension of energy. In fact, this scaling argument 
can be applied to any direction of the applied field. On the other hand, the explicit form of 
the scaling function / cannot be determined by the scaling argument alone. 

Now let us calculate the linewidth explicitly for the anisotropy parallel to the applied 
field. As we have discussed, all the effect of the applied uniform field is represented by the 
shift of the field (3.8). Consequently, the perturbation under the applied field H is 

c a = - ¥L(j> R + j' L ) - ^f. (4.3) 

The third term is a constant and thus can be ignored. The second term is 

XH {r R + Jl) = -2nXH' d * (4.4) 



V2 R V^dx' 

which is equivalent to the additional magnetic field of —2n\H. This can be absorbed by 
a renormalization of the magnetic field, giving the shift of the resonance by —27r\H. This 
shift is first order in the perturbation 8 and the field H. 

Now, let us discuss the effect of the first term. We should calculate the correlation 
function ( J + J~) in the presence of the perturbation —\J Z R J Z L . For this particular problem, 
this can be done exactly, because the perturbation J R J[ is proportional to the kinetic term 
of the free boson Lagrangian; it just gives a renormalization of the compactification radius. 
That is, the Lagrangian density reads 

£ = \(M) 2 - MIJZ 

= (4.5) 

Rescaling the field so that the coefficient of the kinetic term is again given by 1/2, the 
renormalized radius R is given as 
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We note that, we have not included the similar renormalization due to the applied field 
because of the subtleties explained in Section III A. In contrast, the exchange anisotropy 
does break the £77(2) symmetry; there is no reason not to include the renormalization in 
the present case. 

The conformal weight of the vertex operator J ± = e ± i2 «^+ i <t>/R j s (A, A) = (1+A', A') or 
(A', 1+A') where A' = (2ttR—1/R) 2 / (8ir) ~ ir 2 \ 2 . Its Green's function at finite temperature 
is given in eq. (3.19). As explained, the Green's function evaluated at the momenta ±77 is 
relevant for ESR. Near the center of the resonance, the spectrum is dominated by the pole 
of the T function; it reduces to 

/x const. 

Thus the resonance is Lorentzian with the width 

r] = AnA'T = 4tt 3 A 2 T. (4.8) 

This is consistent with the scaling argument (4.2). 

To summarize, the exchange anisotropy with the axis parallel to the applied field gives 
the following effects on paramagnetic ESR. 

shift -2ttAT7 oc -H5 
width 4tt 3 A 2 T oc {5/JfT 



V. SELF-ENERGY APPROACH 

In the last Section, the ESR absorption spectrum was calculated directly in the low- 
energy effective theory. This was made possible because the effective theory reduced to the 
free boson theory. However, in general, the problem is more difficult because the effective 
field theory involves interactions. 

A possible application of the field theory method to ESR is to evaluate the Green's 
function appearing in MK formula (2.7) by means of the field theory. While the MK formula 
has been applied to quantum spin chains by several authors, most of the calculations are 
based on classical or high-temperature approximations which break down at low temperature 
and in low dimensions. Thus it would be worthwhile to evaluate the MK formula using field 
theory to study quantum spin systems at lower temperature and in lower dimensions. On 
the other hand, the crucial assumption of the (single) Lorentzian lineshape is made in using 
the MK formula usually without a rigorous justification. Moreover, the MK formula ignores 
the possible polarization dependence discussed in Sec. II C. Thus, in this section, we develop 
a new, systematic field-theory approach to ESR, which we call the self-energy approach. The 
ESR spectrum is given by the imaginary part of the retarded Green's function of S ± . As 
we have discussed in the last section, it corresponds to the Green's function of the current 
operators in the effective field theory via eq. (3.27). 

We now assume that the perturbation preserves a symmetry which forbids mixing be- 
tween J x and J y , namely (J x J y ) = 0. Then the correlation function of the total spin can 
be decoupled to J x and J y part. 
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(S + (t)S-(0)} = ^Jd Xl J dx 2 (J x R (t,x 1 )e iH ^-^J x R (0,x 2 )) + (Jl(t,x 1 )e- iH ^- x ^J x L (0,x 2 )) 

+(r R (t,x 1 y H ^ + ^r L (o,x 2 )) + (jz(t, Xl )e- iH ( xi+x *JZ(o,x2)) 

+(J X -> J y ) (5.1) 

Since our effective field theory is S'?7(2)-symmetric, we may freely rotate the xyz-axes. 
Thus, instead of calculating correlation functions of J x we can calculate those of J z , with 
perturbations also rotated correspondingly. The same applies to calculation of J y correlations. 

The motivation for us to rotate the xyz-axes is that, J z is expressed as a derivative 
of the boson field as in eq. (3.21). Thus the problem is reduced to the calculation of 
the bosonic correlation function (00). The structure of the bosonic correlation function is 
well established by the standard diagrammatic perturbation theory, and the ESR lineshape 
is related to the boson self-energy as we will show below. On the other hand, when the 
perturbation allows mixing of J x and J y (in the original representation), there seems no 
way to reduce the problem to the (00) correlation function. In such cases, we do not know 
at present how to construct the theory of ESR based on self-energy. Thus, below we restrict 
ourselves to the situation in which J x and J y do not mix, in the discussion of the self-energy 
approach. We remark that there is no apparent difficulty in the application of the MK 
formula even in cases where the perturbation allows mixing of J x and J y . 

As mentioned in Section III B, we restrict ourselves to the case where the perturbation 
does not contain an oscillating factor e lHx . Then the contribution from the cross terms such 
as (JrJl) vanish in eq. (5.1), due to momentum conservation. The correlation function thus 
reduces, upon Fourier transformation to 

(S+S-)(cu) = JL [(J x R J x R )(co, -H) + (J x l J x l )(lu,H) + («)(c, -H) + (JIJI)( U ,H)] (5.2) 

where (JJ)(u,q) denotes the correlation function at frequency u and momentum q. As we 
have discussed above, we now rotate the axes and calculate J z correlation function instead 
of J x and J y , to obtain 

(S+S-)(u,) = 

^L[(J R J R ) X ^, H) + (JlJD^u, -H) + {J R J z R ) y ^ z {u, H) + (JlJl) y _> z (tu, -H)}, (5.3) 

where () x _> z means the correlation function with the perturbation rotated x — > z. Using 
eq. (3.21) and (3.23), those correlation functions can be written in terms of bosonic correla- 
tion function: 

(S+S-)(cu) = t±JpL{M)^,( u ,H) + iLU+ 4 f )2 (H) y ^(co,H), (5.4) 

where we have used the symmetry (00) (uj, —H) = (00) (cu, H). The above formula is useful 
if the perturbation (after the rotation) is given by a Lagrangian density local in the boson 
field 0. If, for example, the Lagrangian density is local in terms of the dual field after the 
rotation y — > z, the second term in eq. (5.4) should be replaced by 

U ' + //)2 (00)^(c,iJ). (5.5) 



47T 
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In fact, there is a subtlety in defining the current. In the free boson theory without 
interactions, we have 



d^ = W (5 - 6) 

dt~ &c' [ } 

so that we may represent the current operator in terms of either or 0. However, in 
the presence of the interaction, we cannot define the dual fields and that satisfy both 
identities. For example, let us take the Lagrangian density 

C = \(d^) 2 -Acos/30, (5.8) 
and define the dual field by eq. (5.6). Then, from the equation of motion, we find 

d x <j>(t, x) + d t (f)(t, x) = -/3\ f cos /3<f>(t', x)dt', (5.9) 

violating eq. (5.7). 

Thus it is not completely clear whether the current operator should be written as a 
derivative of or 0. However, upon Fourier transform, the "difference term" (right-hand side 
of eq. (5.9)) does not give a sharp peak. (Recall that only the operators of conformal weight 
(1, 0) or (0, 1) produce a delta-function spectrum. Other operators give broad spectrum given 
by eq. (3.19), even in the zeroth order.) Moreover, the contribution from the difference term 
is suppressed by a factor A 2 . Therefore, the difference term would lead, at most, only to 
a small and broad background. In discussing the lineshape of the main resonance, we can 
ignore the difference term and focus on the derivative of either boson field or 0. For 
calculational convenience, we choose to use (or 0) if the interaction is given in terms of 
(00 

Thus the problem of finding the ESR absorption spectrum is reduced to the calculation 
of the correlation function of the boson field 0. We now make the Wick rotation and consider 
the corresponding Matsubara Green's function defined by 

Gab (t) = ~\™r [Mr)B(0)} , (5.10) 

where T T is the ordering operator with respect to the imaginary time r, and A(t) = 
e T7i Ae~ T7i . The standard diagrammatic perturbation theory can be applied to the Matsub- 
ara Green's function. After obtaining the Matsubara Green's function, we can analytically 
continue back to real time to obtain the retarded Green's function. 

Provided that the Lagrangian is local in terms of the boson field, its correlation function 
can be written in a self-energy form: 

Q H {uJ n ,q) = — ^— r, (5.11) 

ul + q 2 + U(u n , q) 
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where Q is the (full) Matsubara Green's function, u n is the Matsubara frequency, and 
H(w n , q) is the self-energy, namely the sum of all one-particle irreducible diagrams. Thus we 
obtain 

{iu n + Hf -1 (^ n + #) 2 -1 



47T ^ + # 2 + IL^, #) 4?r cU 2 +# 2 + rL> n ,#)' 

(5.12) 

where ITa, and IT y are the self-energy in the Matsubara formalism, respectively for (<fxf)) x ^ z 
and {(f>(f))y^ z . This gives, upon the analytic continuation, the retarded Green's function 

R (uj + H) 2 1 (io + iJ) 2 1 

4tt c, 2 -tf 2 -II^,tf) + 4tt ^-H^-U^HY (5 - 13) 

where the "self-energy" Ilf (a = x, y) is defined by the analytic continuation 

n2(iu n ,q) = Il a (u) n ,q) (5.14) 

for uj n > 0. 

First let us check what we obtain in the absence of the perturbation. Then = IT^ = 
so that the Green's function has a pole at u — H: 

Gf +(J _(w)~ — k . (5.15) 

This means that we have a completely sharp resonance at the Zeeman energy uj = H as 
expected, in agreement with the equation of motion. The residue H/tt at the pole of the 
Green's function gives the intensity of the resonance. This is also consistent with the exact 
result from the original spin chain. 

„ r°° 9m 



;) = -ij o °°dt([S + (t),S-(0)}) = Q , (5.16) 



where m is the magnetization. For small field H, the magnetization is given by m = XuH, 
where the uniform susceptibility is 

*« = ^ < 5 - 17 > 

in the low-temperature limit, ignoring the effect of the isotropic marginal operator 29 . (We 
remind the reader that we have been setting v — 1.) Thus we obtain the amplitude 2m = 
H/n, in agreement with eq. (5.15). 

A symmetry breaking perturbation 7i' would give no n- vanishing boson self-energy U x , U y . 
This changes the ESR lineshape. Near the resonance uj ~ H, we can write 

H 1 H 1 

Qs+s-M = tou-H-^U^H) + tou>-H-±H*(u>,H) (5 ' 18) 

If the self-energy changes smoothly around the resonance u ~ H, we may regard the self- 
energy as being constant in a frequency range sufficiently close to the center of resonance. 
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Then, within this range, the lineshape is given by a Lorentzian, and the real and imaginary 
parts of the self-energy give the shift and width of the ESR, respectively. The linewidth is 
given by 



-1 



imn*(#,#), 



(5.19) 



7/ 



2H 



while the shift is 



1 



(5.20) 



2H 



for a — x,y. In general, the signal could be superposition of two Lorentzian spectra corre- 
sponding to LT^ and LT^. However, in the concrete cases we study in the present paper, IT^ 
and LT^ are equal; thus a single Lorentzian lineshape is predicted. 

Therefore we have successfully formulated the theory of ESR without any particular 
assumption on the lineshape. The self-energy is usually a smooth function of uj near uj ~ H 
for finite H except for the smooth weak background discussed below Eq. (5.9); we have 
given a microscopic foundation for the Lorentzian lineshape which is assumed a priori in 
the MK approach. Application of the present self-energy formalism to two cases relevant 
to experiments will be discussed in the following sections. However, precisely speaking, our 
approach is only formulated ignoring the isotropic marginal operator, which is generally 
present in the effective theory of the Heisenberg antiferromagnetic chains. Some discussions 
on the effects of the isotropic marginal operator will be given in Section VIC. 

Comparing with the assumption (2.11) used in our derivation of the MK formula in 
Appendix A, it is obvious that the MK formula and the self-energy approach are closely 
related. Namely, S introduced in eq. (2.11) corresponds to n R /(2i7) if they vary smoothly 
around the resonance. The important difference is that it is an assumption that the Green's 
function can be written as in eq. (2.11) with a smooth E whereas we can prove eq. (5.11) 
using the diagrammatic perturbation theory. The self-energy, n, is given by the sum of all 
one-particle irreducible Feynman diagrams as in proven in any book on field theory. In this 
way, our self-energy formulation effectively gives a proof of the Lorentzian form (2.11) which 
is often assumed without a microscopic foundation. We emphasize that, although eq. (2.11) 
may appear innocent, it is a rather strong assumption and is far from trivial. When the 
lineshape turns out to be Lorentzian, the results must agree between the MK and self-energy 
approaches, if the correlation functions are evaluated correctly. This will be verified for a few 
cases in Sections VI B, VIII B and VIII C. On the other hand, while the validity of the MK 
formula is limited to the lowest order perturbation theory, the self-energy formulation allows 
us to go beyond that. In fact, we will make a non-perturbative analysis of the lineshape, 
based on the self-energy formalism, in Section VIII E. 

We note that assumptions similar to eq. (2.11) have been made in the literature for 
different problems. Sometimes the (counterpart of) assumed S in eq. (2.11) is referred to as 
the memory function. For example, Giamarchi 30 studied the conductivity of the TL liquid 
with the bosonization method. His discussion was rather closely related to our analysis of 
ESR in the present paper. (See also Ref. 32.) In fact, he calculated the ac conductivity of a 
TL liquid by evaluating the memory function with the field theory. This is quite similar to a 
field-theory calculation of the MK formula for ESR, which we will discuss in later sections. 
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We could also apply our self-energy approach to problems such as the conductivity of 
TL liquids. This might be useful for providing a more rigorous justification and a possibility 
to go beyond the lowest order perturbation theory. The possible breakdown of the MK 
formula, in the context of the conductivity of a TL liquid, was discussed by Giamarchi and 
Millis 31 . 



VI. EXCHANGE ANISOTROPY PERPENDICULAR TO THE MAGNETIC 

FIELD 

Now we consider the exchange anisotropy with the axis perpendicular to the applied 
magnetic field. Let us take the axis of the anisotropy as the x-axis. In the low-energy 
effective theory, at zero uniform field, the anisotropy term is given as 

r — — \ t x t x 

= -\{J X R J1 - JPD + \JrJI + \j R ■ Jl- (6.1) 

Here the parameter A, which is proportional to 5 for a small 5, is the same as the one intro- 

— # — * 

duced in eq. (4.1). The last term Jr- Jl of the second line is the isotropic marginal operator, 
which does not affect the resonance directly and will thus be ignored in the following. 

Now let us include the effects (3.8) of the applied uniform field H. The first and second 
terms in eq. (6.1) are transformed into 

A = ~{J x rJ1 ~ JVl) + \JrJI + ^(J R + Jl) + ^f. (6.2) 

Fortunately, there is no oscillating factor e lHx here. The last constant term has no effect in 
the following, and will be ignored. The third term represents the additional magnetic field 
of +ti\H. (Compare with eq. (4.3).) This is again absorbed by a renormalization of the 
uniform field H, giving the shift of 7r\H. 

The remaining problem then is to study the effect of the perturbation 

^ = -\wI-JIJI~JrJI) (6.3) 

The first two terms corresponds to an interaction in terms of the boson field 0, and the 
problem cannot be reduced to a free field theory. Thus it is not possible to calculate the 
ESR absorption spectrum directly as we have done for the exchange anisotropy parallel 
to the magnetic field in Section IV. Therefore, we will employ the self-energy approach 
developed in Section V. 



A. Self-energy approach 

Because the anisotropy considered here breaks the rotational symmetry in the xy-plane, 
we expect a polarization dependence. Thus let us consider the correlation function of S x 
and S y separately. Under the magnetic field, S x ' y at zero momentum are expressed as 
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JUH) + Jl(-H) + J R ( -H) + J L (H) ^ 



2V8 



7T 



S y = Jr(H) + J£(-H) - J R (-H) - Jj(H) 

We emphasize here that, under the magnetic field, S x is related to both current operators 
J x and J y . The original spin operator and the current operator are quite different objects. 

Absorbing the third term in (6.2) as a renormalization of the magnetic field, the pertur- 
bation respects the symmetry J x — > J x , J y — > — J y , J z — > —J z . Thus the cross term (J x J y ) 
vanishes in this case, allowing us to proceed with the rotation trick described in Section V. 
Namely, 

(s*s*) = ^[{JrJr) + WD + WD + WD 

+WD + WD + WD + WDKq = h) 

~\{JrJr)x + (JrJDx + (JlJr)x + {JlJDx 



32tt 2 

+ <«>» + WDv + (^^>vl(9 = H) 
+^[..](q = -H), (6.6) 

= IWr) - WD - Wr) + WD 
+WD - WD - WD + WD](q = h) 

+ [...](q = -H) 

= [(JrJDx ~ (JrJDx ~~ (JLJDx + {JlJDx 
+WDv - WDv - WDv + WDvKq = h) 

+ [...](q = -H). (6.7) 

(s x s y ) = ^IWD - WD + WD - WD 

+WD + WD - WD - WDKq = h) 

-\..]{q = -H) 

= [{JrJDx ~ {Jr^Dx + (JL^Dx ~ (Jl-JDz 
+WDv + WDv - WDv - WDvKq = h) 

-[.■■](? = -H). (6.8) 

Here () x and () y means the expectation value in the presence of the (rotated) perturbation 
X/2[J[J Z R - (J X L J X R + J y L J y R )} and A/2[( Jf J X R - J y L J y R ) - J[J Z R }, respectively. Fortunately, these 
can be written in terms of either or 0: 

J'lJ'r ~ Wr + JIJr) = -Ad,<P) 2 ~ 2 cos Vto<f>, (6.9) 

Wr ~ JPD ~ JIJr = M^0) 2 + 2 cos ^0. (6.10) 

The (d^<f)) 2 term gives a renormalization of the radius R. However, in the lowest order of 
the perturbation theory, its effect is negligible on the boson correlation function (00) and 
thus will be dropped in the following. 
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Thus, in evaluating (J Z J Z ) X we will represent the current operator J z as a derivative of 
0, so that the problem is reduced to the correlation function of the fundamental boson field 
(f) in the presence of the interaction in terms of <fi. On the other hand, in evaluating (J z J z ) y , 
we will express the current J z by of for (J z J z ) y . 

As a result, we have 



nR '■> 



7T 



= IT 
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I 7Z IZ 

J R J L 
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(6.11) 
(6.12) 
(6.13) 
(6.14) 
(6.15) 
(6.16) 



where Q R (a = x, y) is the retarded Green's function defined by the expectation value 
(. . .) Q , U r (uj, q) is the self-energy for the boson field <fi in the presence of the interaction 
—A cos \^87T(f) (or the self-energy for the boson field in the presence of A cos y/Sn^, but this 
is identical.) Plugging these into eqs. (6. 6), (6. 7), we obtain 



->R 



H 2 



2tt uo 2 


-H 2 ~U R (u,HY 


uj 2 


1 


27TL0 2 


- H 2 - U R (cu, H) ' 


uoH 
i 


1 



2tt lo 2 - H 2 -U r {uj,HY 



(6.17) 
(6.18) 
(6.19) 



where Q R p is the retarded Green's functions of the spin operators S a and S 13 , as defined in 
eq. (2.3). 

For a direction a in the xy-plane, 



■>r 



H 2 cos 2 $ + lu 2 sin 2 $ 



7T 



u 2 - H 2 -U R {u,Hy 



(6.20) 



where $ is the angle between x and a directions, namely the angle between the anisotropy 
axis and the polarization of the electromagnetic wave. 

As a result, for any directions of the polarization perpendicular to the magnetic field, the 
ESR lineshape is Lorentzian with the width — lmU R (H, H)/(2H). However, the lineshape 
has some angle dependence through the numerator H 2 cos 2 9 + u 2 sin 2 9. In fact, the present 
result is consistent with the exact and rigorous relation (2.18) for original spin model. This 
serves as a consistency check of our field-theory approach. 
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Now let us calculate the self-energy II of boson field in the presence of interaction 
A cos \/&K(j). It is easy to see the first order perturbation to the boson correlation function 
vanishes due to symmetry. The second order perturbation to the boson correlation function 
does not vanish and can be calculated by the diagrammatic expansion (ie. Wick's theorem). 
The second-order term in the boson correlation function is related to 

* m y^) e -^m m = ^J2 ( ^ rHv/ f r (0(l) = 0"(2) :: <T (3) : 0(4)) (6.21) 

= ttA 2 (0(1)0(2)) ( e iV5^(2) e -iV§^(3)) (0(3)0(4)) + (2^3) 
-7rA 2 (0(l)0(2))(e i ^ (2) e- l ^ (3) ) (0(2)0(4)) + (2^3) 
+ttA 2 (0(1)0(4)) ( e i^(2) e -i^(3)j ( 6>22 ) 

The three terms here represent contributions form different kinds of Feynman diagrams, as 
shown in Fig. 3. The second type of the term -(0(1)0(2)) ^ e W^4>m e -W^m) (0(2)0(4)) + 
(2 <-> 3) represents the "tadpole" type Feynman diagram (Fig. 3 (b)), while the last term 
corresponds to a disconnected Feynman diagram (Fig. 3 (c)), which is canceled by the 
correction to the partition function. 



(a) 



(b) 




(c) 




FIG. 3. Three types of Feynman diagrams appearing in the perturbative expansion. (a),(b) 
and (c) correspond to the first, second and third terms in eq. (6.22), respectively. The disconnected 
diagram (c) is canceled by the correction to the partition function; the "tadpole" diagram (b) does 
not contribute to the imaginary part of the self-energy (i.e. the linewidth) 
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In fact, there is a similar contribution from e -*V8n<P(2) e +iVsn<t>{3) fogging the above, and 
one has to integrate the coordinates 2 and 3 over Euclidean space-time. As a result, we 
obtain the self-energy in the lowest order (0(A 2 )) of the perturbation as 

U(u n , q) = 47rA 2 [G (lil) K, q) - G ( i,i)(0, 0)], (6.23) 

where C(i,i) is the Matsubara Green's function of the operator of the conformal weight (1, 1) 
in the free boson theory. These two terms come from type (a) and (b) Feynman diagrams 
in Fig. 3, respectively. Analytic continuation back to real time leads to 

= 47rA 2 [Gf lil) (^,g) - 0^(0,0)], (6.24) 

where is the retarded Green's function corresponding to the Matsubara Green's func- 

tion G(i,i). Its imaginary part can be derived by taking the limit A, A — > 1 in eq. (3.19): 



lm[-G? hl) (uj,q)] = ?-(u J 2 -q 2 ) 



coth 10 ' + coth — 

4T 4T 



(6.25) 



8 

The imaginary part then reads 

-ImU R (H, H) = 47r 3 X 2 HT, (6.26) 

giving the width 

7] = 2tt 3 A 2 T (6.27) 

Again, this is consistent with the scaling analysis (4.2). The real part is proportional to 
{u 2 — q 2 ), which corresponds to a wavefunction renormalization, and does not lead to any 
shift at 0(A 2 ). In any case, there is a shift of 0(A) discussed above, which is dominant. 

To summarize, the exchange anisotropy with the axis perpendicular to the applied field 
gives the following effects on paramagnetic ESR. 

shift +nXHocH5 

width 2tt 3 A 2 T oc {S/J) 2 T 

Comparing to the result for the exchange anisotropy with the axis parallel to the applied 
field, the width obtained here is half of the result (4.8) for the parallel case. This can be 
understood naturally with the MK formula as we will discuss in the next subsection. On 
the other hand, the shift takes opposite sign and the absolute value is half of that in the 
parallel case. 



B. MK approach 

The lineshape is shown to be Lorentzian in the two cases discussed above (exchange 
anisotropy parallel and perpendicular to the applied field), up to a possible broad background 
of 0(A 2 ). Thus the MK formula is expected to be also valid for these cases. In order to 
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check consistency of our field-theory approach, here we study the same problem with the 
MK formula. 

Let us consider the exchange anisotropy parallel to the applied field considered in Sec- 
tion IV. We may apply the MK formula to the spin chain Hamiltonian and then take the 
continuum limit, but taking the continuum limit first and then apply the MK formula turns 
out to be simpler. Absorbing the second term of the effective perturbation (4.3) into a 
renormalization of the magnetic field, we need to consider the effect of the perturbation 
W ~ A / dxJ[J R . 

First we have to obtain the commutator (2.9) appearing in the MK formula. The total 
spin raising/lowering operator S ± in the continuum limit is given from eq. (3.27) as 

S± = -jL= J dx(J±e ±iHx + Jte***). (6.28) 

Using the standard commutation relation among the currents, the commutator A is given 
by 

A = i[H', S + ] =i\J dx[Jl{x)J+{x)e lHx + J z R (x) J+(x)e- iHx ]. (6.29) 

J[J R and J R J£ are primary fields with the conformal weight (1, 1). Thus, from the MK 
formula (2.7) we obtain the linewidth 



2A 2 



v =^-Im[-G? hl) (H,H)} (6.30) 



XuH 



r (i,i) 



The Green's function is what we have already considered in (6.25), and thus we obtain the 
width 

r] = — tt 2 T. (6.31) 

Xu 

Using eq. (5.17) again (recall we have set v = 1), 

7] = 4tt 3 A 2 T. (6.32) 

This indeed agrees exactly with the result (4.8) obtained by quite a different approach. We 
remark that a similar derivation of a similar formula for the ac conductivity of a TL liquid 
was given earlier by Giamarchi 30 . 

Next let us consider the exchange anisotropy perpendicular to the applied field. Absorb- 
ing the third term in (6.2) into the renormalization of the magnetic field, the perturbation 
to be considered is H! = (A/2) / dx(JlJ R — J V L J V R — JIJr)- Consequently, the commutator 
becomes 



A = i[H', S + ] =i±J dx[(J£(x)J R (x) + J z L JtY Hx + (Jl(x)J R (x) + JtJ z R )e- iH % (6.33) 



.A 

This leads to 
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rj = i m [_ G « i} {H) H)] = 2vr 3 A 2 T, (6.34) 

where we have used the susceptibility (5.17) in the second equality. Again we have found an 
exact agreement with the self-energy approach (6.27). The ratio 2 of the width between the 
parallel case (4.8) and the perpendicular case (6.27) is simply understood in this approach. 
It arises from the factor of 1/2 and the presence of twice as many terms in Eq. (6.33) as 
compared to Eq. (6.29). In fact, such an angle dependence also holds at higher temperature 
and has been discussed in the literature, for example in Refs. 33, 34. 



C. Effect of the marginal isotropic operator: logarithmic correction 

The Hamiltonian of the Heisenberg antiferromagnetic chain with a small anisotropy in 
the z direction can be written as 

n = n - [ g x (r R r L + jpi) + g z r R r L ] , (6.35) 

where we ignored the applied field H, which will be considered later. Here we can rewrite 
the perturbation as 

g x J L -J R + (g z -g x )J z L J z R , (6-36) 

where the first term is the isotropic marginal operator. The second term gives the anisotropic 
interaction A = — g z + g x . 

As is now well known, the isotropic marginal perturbation exists in the low-energy ef- 
fective theory of the Heisenberg antiferromagnetic chain, giving several effects such as the 
logarithmic correction to the magnetic susceptibility 29 at low temperature. While it has a 
simple form J L • J R at H — 0, it becomes complicated if we include the effect of the ap- 
plied field H. It introduces complications such as the momentum non- conservation in the 
effective theory and the mixing of J x and J y , thereby invalidating the simple self-energy 
approach discussed in Section V. Thus we actually have no microscopic derivation of the 
Lorentzian lineshape in the presence of the isotropic marginal operator, at present. On the 
other hand, the operator by itself, being isotropic, does not directly affect the linewidth. 
Since the isotropic marginal coupling constant renormalizes to zero, we may expect the 
Lorentzian lineshape is basically unaffected by its presence. It does, however, indirectly af- 
fect the linewidth through the renormalization of the anisotropic perturbation as we discuss 
in the following. 

As discussed in Ref. 42, the coupling constants g x and g z are renormalized by the 
Kosterlitz-Thouless type RG flow. The solution of the RG equation (for H = 0) in the 
lowest order gives 

f = T- • w 1 , \ (6-37) 
47rsmh(elnr) v ' 

g z = -^coth(elnr), (6.38) 

4"7T 
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where r is the scale variable (oc J/T) and e is a constant, which determines the crossover 
scale. [This solution is valid only if the infrared (IR) limit is a massless free boson theory, 
namely if 5 < 0. We proceed by assuming this case; the final result on the ESR linewidth 
should be valid also for S > 0.] In the IR limit r — > oo, g x = and 

«f (oo) = ± (6.39) 

This corresponds to a renormalized free boson Lagrangian (1 — 2irg z (oo))(d fl (j)) 2 /2, which 
leads to the critical exponent r/ z = 1 — 2irg z (oo), where (S z (r)S z (0)) ~ r~ Vz . 

On the other hand, the critical exponent in the low-energy limit of the Heisenberg XXZ 
model has been obtained from the Bethe Ansatz exact solution. For the Heisenberg model 
with an exchange anisotropy 

H = £ J{S%S% +X + S]S] +1 ) + ( J + 6)S Z S Z +1 , (6.40) 

3 

it is known that 

% = ^s = l--cos- 1 [l + 4 (6-41) 
for a negative 5. Combining these results, we obtain, for small e, 5, 



-85 



J 



(6.42) 



Since the isotropic part g x J2 a JlJr commutes with S + , the important perturbation is the 
"asymmetric part" g z — g x ■ In the intermediate scale r e 1 ^, which would be relevant to 
ESR for a weak anisotropy, 

X = -g z + g x = ^e 2 lnr = (6.43) 

87T 7T d J 

This corresponds to the coefficient A introduced in eq. (4.1). The larger of the temperature 
T or the applied field H imposes the cutoff of the RG flow, and thus the scale factor r should 
be replaced by J/max(T, H) . 

In the present discussion, the uniform field H appears only as a cutoff scale imposed on 
the RG flow at zero field. Thus, to this order, the renormalization of the coupling constant 
A applies to arbitrary direction of the anisotropy relative to the applied field. Therefore we 
conclude the low-temperature asymptotic behavior of the linewidth and shift to be 

^ = -v^^Tjry (6 - 45) 

if the anisotropy axis is parallel to the applied field. They are 
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j V 



max(T, H) 



T, (6.46) 



^ = - 2 -Jn J —— y (6.47) 

n 2 J max(T, H) 

if the anisotropy axis is perpendicular to the applied field. The shift depends on the sign 
of the anisotropy. When comparing with experiments or existing literature, it should be 
recalled that we discuss the shift in frequency (for a fixed field H) while usually a shift 
in the resonance field for a fixed frequency uo is studied. For example, in the presence of 
the dipolar interaction, which corresponds to negative 5, when the field is applied parallel 
to the chain (ie. anisotropy) axis we obtain a positive frequency shift, namely a negative 
shift in the resonance field. The shift is in the opposite direction when the applied field 
is perpendicular to the chain axis. These conclusions are qualitatively consistent with the 
literature. 48,14 



D. Comparison with experiments 



In this paper, we have not calculated the ESR lineshape for a general relative direction 
between the anisotropy axis and the magnetic field, let alone more complicated anisotropy 
of general form. However, the results (4.8), (6.27) together with the scaling argument (4.2) 
imply that the linewidth due to the exchange anisotropy (or dipolar interaction) scales 
proportionally to the temperature T in the low temperature regime T <^ J (but above 
the Neel or spin-Peierls transition temperature). This, in fact, appears to be observed in 
many quasi-one dimensional S = 1/2 antiferromagnets 35-37 ' 34 ' 39 ' 40 including CPC, KCUF3, 
CuGeC>3 and NaX^Os. In the case of Cu benzoate 10 , there is a field-dependent diverging 
contribution to the linewidth at low temperature due to a staggered field effect, as we will 
discuss in Section VIII. There seems to be another contribution to the linewidth, which is 
approximately T-linear and frequency-independent. We presume the latter contribution is 
due to the exchange anisotropy. 
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FIG. 4. The temperature dependence of the ESR linewidth in KCUF3, CuGe0 3 and NaV 2 0s. 
The data are taken from Refs. 34, 39, 40, respectively. The horizontal axis is the temperature T 
normalized by the exchange coupling J, and the vertical axis is the normalized linewidth. 



In Fig. 4 we show the observed 34 ' 39 ' 40 ESR linewidth for KCuF 3 , CuGe0 3 and NaV 2 5 , 
as a function of the normalized temperature T/J. We note that, these materials exhibit 
phase transitions (such as Neel and spin-Peierls transitions) at low enough temperatures, 
where the linewidth appears to diverge. Since we focus on one-dimensional systems in the 
present work, in Fig. 4 we have omitted such temperature regimes, above which we may 
regard the system simply as a spin chain. It could be possible that, however, the displayed 
data are still affected by the interchain interactions, the spin-Peierls instability etc. 

An analysis on the linewidth in NaV 2 05 similar to ours was published previously by 
Zvyagin 41 . However, we also remark that the T-linear behavior of the linewidth due to an 
exchange anisotropy was reported earlier in Ref. 12. In fact, eq. (3) in Ref. 41 is equivalent to 
eq. (11) in Ref. 12. Moreover, in Ref. 41 it was argued that a bond-alternation perturbation 
leads to a linewidth oc ( J 1 — J 2 ) 2 /T 2 . However, the argument (leading to Eq. (4) in Ref. 41) 
cannot be correct per se, because the ESR linewidth must remain strictly zero as long as all 
terms in the Hamiltonian except the Zeeman term commute with the total spin operators, as 
we reviewed in Sec. II A. An isotropic bond-alternation has this property. It is possible that 
an isotropic bond-alternation perturbation (J x — J 2 ) together with an anisotropic uniform 
exchange perturbation 5 might lead to a width, but it would be suppressed by the factor 
S 2 /J 2 as the width should vanish when 5 = 0. In any case, a reliable derivation seems 
lacking so far. We point out that the ESR spectrum cannot simply be related to the boson 
propagator in the field theory, in the presence of a bond-alternation. (See remarks below 
eq. (5.1).) 

In Fig. 4 we took J = 400K, 150K and 560K respectively 34 ' 39 ' 40 for KCuF 3 , CuGe0 3 and 
NaV 2 5 , while there are some uncertainties in the estimate. The linewidth is renormalized 
to be compared with T/J. The low-temperature asymptotic behavior of the linewidth 
indeed seems consistent, although not perfectly, with the universal T-linear behavior we have 
derived. On the other hand, it is difficult to discuss the predicted logarithmic correction in 
the present data. Regarding Fig. 4 as a fitting, the low-temperature asymptotic behavior 
reads 



— ~ < 
T 



' 4.2 x 1(T 4 (CuGe0 3 ,# || c) 
4.7 x 1(T 4 (CuGe0 3 ,# || a) 
17xlCT 4 (KCnFv.HWc) 



17xl0- 4 (KCuF 3 ,tf||c) 

22xl0~ 4 (KCuF 3 , H \\ a) ' { } 



1.3 x 10" 4 (NaV 2 5 ,# || c) 
0.65 x 1(T 4 (NaV 2 5 ,# || b) 



which are given as dimensionless numbers. In these materials, the data for H || a and H || b 
are quite similar, and thus only one set of them is shown for each material. 

Comparing with our results (6.44) and (6.46), the anisotropy 6/ J seems to be about a few 
percent. It was argued 34 ' 39 ' 40 that, in these material it is too (up to 10 times) big compared 



33 



to what we expect from Moriya's 27 estimate 5 ~ (Ag/g) 2 J where Ag is the anisotropy of 
the g-tensor. (Actually the discussion in Refs. 34, 38-40 was based on the high-temperature 
limit. See Section IX for relation to our low-temperature theory.) However, we believe that 
Moriya's formula is only valid as an order-of-magnitude estimate. There is a room for a 
factor which is presumably not too much different from 1, but could still allow the exchange 
anisotropy that is consistent with the observed linewidth. 

The linewidth deviates from the field theory result oc T at higher temperatures. This is 
not surprising, since the field theory is only valid in the low temperature T <^ J. We will 
give more discussion on the crossover to the high-temperature regime in Section IX. On the 
other hand, if all the materials can be regarded as standard Heisenberg antiferromagnetic 
chains with the same type of anisotropy, we would expect the linewidth to be a universal 
function of T/J. However, in Fig. 4 it is evident that the linewidth behaves differently at 
high temperature, especially in KCuF 3 . This suggests that not all of them can be described 
by the standard Heisenberg antiferromagnetic chain (3.1) with the same type of anisotropy. 
We remark that the low-temperature asymptotic behavior should be universal for a certain 
class of Hamiltonians, but the explicit coefficients obtained in eqs. (6.44) and (6.46) are 
specific to the standard Hamiltonian (3.1). 

Certainly, there are many questions still to be understood. An important problem is 
the dependence on the direction of the applied field. In the case of NaV 2 5 , the observed 
linewidth at low temperature is twice as large when H \\ c compared as when H _L c. This is 
consistent with our result, if an exchange anisotropy with the single anisotropy axis parallel 
to c is assumed. However, in the case of CuGeOs and KCUF3, the observed linewidth for 
H || c is smaller than that for H || a and H || b. This kind of angular dependence cannot be 
explained with an exchange anisotropy with a single anisotropy axis. This suggests that we 
have to consider more general types of anisotropy, or some other effects. 

A complete theoretical description of the experimental data in these materials is left 
for the future. Nevertheless, we believe that the universal decrease of ESR linewidth at 
low temperatures in S = 1/2 antiferromagnetic chains is basically understood with our 
theory. Ours is presumably the first 12 microscopic derivation of this approximately T-linear 
linewidth. In Refs. 34, 38-40 a completely different interpretation was proposed. However, 
we will argue against it in Section IX. 

VII. ESR IN AN XXZ ANTIFERROMAGNET 

So far in this paper, we have restricted ourselves to the case of small anisotropy. However, 
in principle ESR can be measured in a system which is far from isotropic. To apply the 
self-energy formalism to a not small anisotropy, one has to sum up higher orders of the 
perturbation. In addition, the foundation of our self-energy formalism based on the weakly 
broken SU (2) symmetry may be questionable in such cases, because the SU (2) symmetry 
is strongly broken in the spin Hamiltonian. 

However, there is one case in which we can study ESR with a strong anisotropy: an 
easy-plane XXZ antiferromagnet with a field applied perpendicular to the easy plane. This 
is nothing but the isotropic Heisenberg antiferromagnet with a negative exchange anisotropy 
parallel to the applied field (6.40), with 5 < 0. Here we can apply the direct calculation 
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introduced in Sec. IV. 

The compactification radius R for the XXZ model with a given anisotropy 5 is known 
from Bethe Ansatz exact solution and is given in eq. (6.41). Using this radius, the ESR 
absorption spectrum given by the Green's function (3.19) of the vertex operator with the 
conformal weight (3. 17), (3. 18). Since 5 is not small, the spectrum is no longer a simple 
Lorentzian, except at low enough temperature T <^ H <^ J where the spectrum reduces to 
the Lorentzian (4.7). 

In this Lorentzian case, the width here does not reduce to the previous one (6.44) which 
was proportional to 5 2 , even in the limit 5^0. The reason of this disagreement is that they 
describe different regimes. The result (6.44) is valid when the energy scale max(T, H) is 
above the crossover energy E c = e -1 / 6 , while the present result is valid if the relevant energy 
scale T and H are both below E c . For a small anisotropy, the crossover scale is exponentially 
small, making eq. (6.44) realistic for the experimentally accessible regime. 

For a small exchange anisotropy and above the crossover energy E c , the width is pro- 
portional to A 2 in the leading order of perturbation theory; the width is insensitive to the 
sign of the anisotropy (easy-plane or easy-axis). However, when the anisotropy is large or 
T,H <C E c , this symmetry no longer holds. In fact, the system in the zero temperature limit 
is gapless for an easy-plane anisotropy (5 < 0) while it acquires a gap ~ E c for an easy-axis 
anisotropy (5 > 0). In the gapful case 5 > and T,H <C E c , ESR probes the creation of 
the elementary excitation above the groundstate; the absorption spectrum then has a sharp 
peak centered at the energy of the gap. 



VIII. TRANSVERSE STAGGERED FIELD 



As we have discussed in Section III B, a staggered field is the most relevant perturbation 
of the isotropic Heisenberg antiferro magnet. Breaking the SU(2) symmetry, the staggered 
field affects also the ESR spectrum. Here we discuss the effect by the field theory methods 
described in previous sections, and then explain the mysterious observations in ESR exper- 
iments 10,11 on Cu Benzoate in the 1970's which were recently confirmed and extended 45 . 

Let us focus on the case of a transverse staggered field 

w = hj2(-iys*. (8.i) 

j 

As we have discussed already, the staggered field is mapped to the operator 

n x ~ k cos (2rrR4>) (8.2) 

which has scaling dimension 1/2. A standard scaling analysis similar to that in Section IV 
shows that, ignoring the logarithmic correction, the linewidth should be given as 

V = Tg(^,j), (8.3) 

where E g is the excitation gap 6,7 due to the staggered field proportional to /i 2 / 3 J 1 / 3 . Again, 
the scaling argument alone cannot determine the actual form of the scaling function g. 
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A. Self-energy approach 



As we have discussed, The staggered transverse field (8.1) is mapped to the field theory- 
operator n x : 

H' = hJ2(-l) j S* ~khj n x (r)dr (8.4) 
j 

where k is a constant, and we normalize n x by (n x (r)n x (0)) = 1/r. Namely, k 2 gives the 
correlation amplitude (SqS x ) ~ (— \yk 2 / j. This form is not affected by the application of 
the magnetic field H, except for the possible renormalization of the amplitude k and the 
exponent, which we will ignore. 

The SU(2) WZW field theory with the perturbation n x has rotational symmetry about 
the x-axis. While the original spin problem is not invariant under a rotation about the 
x-axis due to the applied field, the effective field theory does have this symmetry. As a 
consequence, correlation functions of the type (J x J y ) vanishes. Thus we can apply the self- 
energy method by reducing the ESR spectrum to Green's function of the bosonic field, as 
discussed in Section V. 

The transverse staggered field in the x direction breaks the rotational symmetry in the 
xy-plane, leading to polarization dependence. Calculations similar to those in Section VI 
lead to the same result (6. 17), (6. 18) and (6.19). The polarization dependence is again 
consistent with the rigorous relation (2.18) which can be applied to the present case. 

The self-energy II is now replaced by the boson self-energy in the presence of the per- 
turbation kh cos \[2t\<\). Again, arguments similar to those in Section VI can be applied to 
obtain the result 

nV,g) = 2n(kh) 2 [Gf 1/4;m (u;,q) - Gf 1/4 , 1/4) (0, 0)], (8.5) 

where the second term comes from the tadpole term. 

The self-energy is a smooth function of u> near the resonance oo ~ H . Thus, the lineshape 
is Lorentzian near the center of the resonance, with the width and shift determined by the 
self-energy at (ou,q) = (H,H). The imaginary part of the second, tadpole term vanishes 
according to eq. (3.19). Using eq. (5.19), the linewidth is given by 

V = — — Im \-G? 1/Atm (H,H) . (8.6) 



H 



From eq. (3.19), the linewidth shows quite a nontrivial dependence on the applied field 
H and temperature T. However, in the weak field regime H <C T, the formula can be 
simplified and linewidth has simple T~ 2 dependence on the temperature. 




1 = T U ™- ( 8 - 7 ) 



We note that this is consistent with the scaling analysis. (Recall that we have set v — 1.) 

The correlation amplitude k 2 was recently determined exactly for the S — 1/2 Heisen- 
berg antiferromagnet 43 ' 42 ' 44 with a logarithmic correction due to the presence of marginal 
operators: 
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The logarithmic correction is translated into a In (J/T) factor in the ESR, where the tem- 
perature gives the IR cutoff. Thus we obtain (upon reinstating v = nJ/2) 

m) ' £ ,„ d) „ 0.685701^ In d). (8.9) 
' 16 V 2 Vr(|) J T 2 y T' T 2 y T> y ' 

Implication of this result on the experiments will be discussed in Sec. V111D. 



B. MK approach 

Since the lineshape is Lorentzian, the MK formula should be valid also in this case, 
provided the correlation function is evaluated appropriately. There are two ways to evaluate 
the commutator (2.9) appearing in the MK formula: to take the continuum limit before 
calculating the commutator, or to first calculate the commutator in terms of the original 
spin variable and then take the continuum limit. We think the former is generally more 
reliable, since the field theory only deals with universal low-energy phenomena while the 
Lorentzian assumption of MK formula would be valid at best in the long-time limit. In the 
present case, the two methods give the same result as we will show below. 

Taking the continuum limit first, we calculate the commutator between the field theory 
operators (8.4) and (3.27). The standard relation between the commutator and OPE leads 
to 

A = [W, S~] 

= [hk J n x (r)dr, JL J J^e'^ + J L (r)e lHr ] 

= (hk/2) J [ e -^<t> e -iHr + e iV2^4> e iHr^ dr ^ 

On the other hand, in the original spin representation, the commutator is easily evaluated 

as 

a = [w , s-] = hj2(-iy s] , (8.ii) 

j 

namely the longitudinal staggered field with the coefficient h. Taking continuum limit, it 
agrees with (8.10). 

Thus, from the MK formula, the linewidth is given by 

V = ^~tM h^/4,1/4) • (8-12) 



Again using eq. (5.17), this agrees exactly with the result (8.6) obtained in the self-energy 
approach. 
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C. Shift of the resonance frequency 



In the present case, there is no shift to first order in h. In fact, the first term in the MK 
formula (2.10) vanishes in the present case. The lowest order shift is thus second order in 
h. This is given by either the MK formula (2.10) or by ReIl R (H, H)/(2H) in the self-energy 
approach. Again, both approaches give the same result for the frequency shift: 

Au = v = -^-Re [-Gf 1/4jl/4) (0, 0) + Gf 1/4>1/4) (1T, H)] . (8.13) 

This is a straightforward consequence of the self-energy approach. On the other hand, 
the derivation from the MK approach might need an explanation. While the second term 
proportional to —G%, 41 i 4 JH,H) just comes from m the MK formula (2.10), the first 

term (proportional to — Gfi/± i/4)(0, 0)) is less obvious. From eq. (8.11), the commutator in 
the first term of the MK formula (2.10) is given by 

[A,S-]=h^2(-iyS- (8.14) 

j 

Its expectation value vanishes if evaluated in the absence of the staggered field H' . However, 
taking the staggered field perturbation into account, 

([A,S-]) = -h 2 Xs + 0(h 3 ), (8.15) 

where Xs is the (transverse) staggered susceptibility. By the linear response theory, we have 

Xs = -k 2 ReGf mm (0,0), (8.16) 

which leads to (8.13), with the replacement of Xu by its zero-temperature limit (5.17). 

For the standard S = 1/2 Heisenberg antiferromagnetic chain, we can apply the exact 
result on the correlation amplitude as we did for the width. We obtain 



1 ' -)V J r, ni) K( ,/n--'^r 



Aoo = -W-ln(-) 



8V2 y T J HT \r(f) 



(8.17) 



For small field H compared to temperature T, we obtain by Taylor expansion of the Gamma 
function 

Th 2 H T 

Au = .344057^^ In (8.18) 

(An incorrect prefactor was given in Ref. 12.) Namely, we obtain the positive shift which 
rapidly increases with decreasing temperature. 

The shift in the presence of the staggered g-tensor was previously discussed by Nagata 46 
using the formula 

Au = -^([[S\H'],S-]). (8.19) 
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derived in Refs. 47, 48. (See also Appendix.) In the present case, it is reduced to the 
expectation value of the staggered field h^^—iyS -. The leading order of the shift in the 
perturbation h is thus given by 



Auj 



(8.20) 



where Xs is the staggered susceptibility 46 . The positive frequency shift (ie. negative field 
shift) was argued to be consistent with the experiment 49,46 . On the other hand, the the- 
oretical result in Ref. 46 is not in quantitative agreement, partly due to the evaluation of 
Xs in the high-temperature classical limit. However, we believe that eq. (8.20) itself is not 
quite correct even if Xs were evaluated exactly. In fact, eq. (8.20) differs from ours (8.13). 
Interestingly, eq. (8.20) is equivalent to including only the tadpole contribution in the self- 
energy approach. The discrepancy becomes particularly important at low magnetic field. 
In the limit of H — > at fixed h and T (although this limit is not realistic in experiment) 
the MK/self-energy approach predict the shift linear in H but (8.20) gives a diverging shift 
~ 1/H, which is presumably unphysical. While eq. (8.19) captures some physics of the 
frequency shift, it fails to include more subtle effects of fluctuation, presumably because of 
the oversimplified ansatz and of not including the long-timescale dynamics. 



The result (8.9) of the perturbation theory implies an interesting behavior of the ESR 
linewidth in materials such as Cu benzoate. As we have discussed in Section III B 1, there an 
effective transverse staggered field is induced proportionally to the applied field {h = cH), 
and the proportionality constant c depends strongly on the direction of the applied field. 
Thus, the linewidth increases as oc T~ 2 as the temperature is lowered. Furthermore, it de- 
pends on the applied field (or the resonance frequency) H and on the direction of the applied 
field. This very characteristic behavior is not expected for the exchange anisotropy. In fact, 
these features were actually observed 10 nearly 30 years ago in ESR on Cu benzoate and ap- 
parently have not been understood until recently. Our results give a natural understanding 
of these observations 12 . 

The only unknown parameters in Cu benzoate were two components of DM vector. We 
have chosen 12 



which seemed most reasonable to fit ESR data 10 . It is also roughly consistent with other 
experiments 7 such as neutron scattering, although not perfectly. This choice of DM vector fit 
rather nicely the direction dependence (Fig. 1 of Ref. 12), temperature and field dependence 
(Fig. 2 of Ref. 12). However, we should note that the determination of the logarithmic 
correction in a practical fitting is a difficult problem; the leading log correction is only valid 
in the low temperature limit. Our fittings were done setting the logarithmic factor to unity. 



D. Comparison with the experiments 



(D a ,i, D c n) = (0.13, 0.02) J 



(8.21) 
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ESR linewidth for H//c" 
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FIG. 5. The temperature and frequency dependence of the ESR linewidth for H \\ c" 10, after 
subtracting the frequency independent part. It is compared our theory (8.9). See also Fig. 2 of 
Ref. 12 for the same comparison without the subtraction. 



There is some discrepancy between the theory and the experiments. We see, in the experi- 
mental data, a field(frequency)-independent contribution which appears to be approximately 
linear in temperature. This is presumably due to effects other than the staggered field. A 
probable mechanism is the effect of an exchange anisotropy, which gives a linewidth which 
is T-linear and independent of the field. If we subtract the field-independent contribution 
from the experimental data at the price of introducing additional fitting parameters, the 
agreement becomes better as shown in Fig. 5. 

Recently Asano et al. made a detailed experimental study 45 on ESR in Cu benzoate. 
They also confirmed our prediction on the linewidth at higher field. In addition, they 
found that, when the temperature is small compared to J but not too low, the shift is 
consistent with our prediction (8.18). See Fig. 2 of Ref. 45. Moreover, we can read off the 
proportionality constant from their Fig. 2 as 

Ato ~ 0.053(— ) 3 (8.22) 

for H || c, where Auo and H are measured in Tesla while T is in Kelvin. On the other hand, 
using the DM vector (8.21), e find that h = 0.095# for H \\ c. Combining this with (8.18), 
we the theoretical prediction 

Auj ~ 0.042(^-) 3 , (8.23) 

where we have again replaced the logarithm haJ/T with unity. Considering the subtlety of 
the logarithmic correction, the agreement between the theory and the experiment is rather 
good. 
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Thus our perturbative results agrees well with the experiments. However, at very low 
temperatures, the lineshape evolves differently than what we expect from the lowest order 
perturbation theory. This will be discussed in the next subsection. 



E. Resonance at very low temperature 

So far, our analysis was perturbative in the staggered field h. While the perturbation 
theory seems reasonable for a small staggered field, it eventually fails at lower temperature 
where the effect of the staggered field is enhanced. In fact, the perturbative expansion turns 
out to be an expansion in Jh 2 /T 3 , which is divergent at low enough temperature. 

The effective field theory describing the S — 1/2 Heisenberg antiferromagnetic chain 
with a staggered field is given by (3.2) perturbed with (8.2). As discussed in Ref. 6, this 
is nothing but the sine-Gordon field theory, which is one of the best understood strongly 
interacting field theories. Since the interaction term (8.2) is relevant, the sine-Gordon field 
theory is massive, ie. has a finite excitation gap E g above the groundstate. The elementary 
excitations of the sine-Gordon model consist of solitons, antisolitons and breathers which 
are boundstates of a soliton and an antisoliton. 

The perturbation theory is expected to be valid only for T ^> E g . Here we consider 
the opposite limit T <C E g , where the system is essentially in the groundstate. Then we 
obtain quite a different picture. It is still valid that the ESR spectrum is given by the 
(00) Green's function at frequency and momentum H. However, we have to consider the 
zero-temperature Green's function in a non-perturbative way. 

In the present case (/3 = V2tt), the lowest excitations are 1st breather, soliton and 
antisoliton, which form an SU(2) triplet. Thus the excitation gap E g is identical to the first 
breather mass M x . The boson field couples to the 1st breather, and thus its propagator is 
given by 



(M)(u,q)~ - — 2 . (8.24) 



Z<S> 

UJ 2_^r_ Mi 

where Z^ is the wavefunction renormalization constant obtained exactly 50 as 



■KV 



Z 4 ' = (l + is) ^exp / dt), (8.25) 

v ; 2sin^ P V Wo sint 



where v = (3 2 /(8-k - /3 2 ). For the present case, i/ = 1/3 and thus Z^ = 0.978689. 

From this, we immediately find that the ESR spectrum at zero temperature is given by 
a delta-function 

(H + y/H* + MA 2 

-ImGg+c-H « A 7 5(cu - JH 2 + M x 2 ). (8.26) 

2y/H 2 + Ml 

Since the wavefunction renormalization Z^ is close to unity, the intensity is identical to that 
of a free resonance. 
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Thus we obtain a rather complicated behavior of ESR in the presence of the staggered 
field. As the temperature is lowered, the linewidth increases in the perturbative regime 
(T ^> E g ) as we discussed, but at lower temperature (T <C E g ) we see a revival of a sharp 
resonance. The width of the resonance vanishes at zero temperature. At small but finite 
temperature < T <C E g , the resonance may be broadened due to the thermally activated 
excitations, but presumably the effect is only of the order of the density of such excitations 



On the other hand, the ESR frequency at zero (or very low) temperature does receive a 
shift due to the staggered field. Namely, the resonance frequency cu is given by 



compared to the Zeeman frequency H. For small mass Mi <C H, the shift is given as 



where we used the result of the breather mass (field-induced gap) E g = Mi ~ 
1.77695[ln(J//i)] 1/6 (J/i 2 ) 1/3 in Refs. 6, 7. Here we emphasize that our self-energy approach 
is valid beyond the lowest order of perturbation theory, unlike the MK formula. Thus it 
allows us a non-perturbative analysis such as the above. 

Our prediction agrees quite well with the experimental result in Ref. 11, as discussed in 
Ref. 12. The non-trivial evolution of the lineshape was indeed observed in the experiment 
in 1970s. Moreover, using the same parameter (8.21) we have used for the perturbative 
analysis, we are able to reproduce the direction dependence of the resonance frequency at 
very low temperature quite well with eq. (8.27) as shown in Fig. 3 of Ref. 12. However, 
the data were only shown at fixed temperature and fixed frequency in Ref. 11. Thus several 
other predictions of our theory could not be compared. After our proposal 12 , Asano et al. 
studied 45 ESR in Cu benzoate at low temperature and at higher field. They confirmed the 
crossover to the non-perturbative regime, and that the resonance at very low temperature 
agreed with the prediction (8.28) for various fields. Moreover, the crossover between the 
perturbative and the non-perturbative regime occurs at temperature T ~ E g , consistently 
with our picture. The broadening at the non-perturbative regime was also consistent with 
the picture e~ Ea ^ T 

On the other hand, the precise lineshape at the crossover temperature regime T ~ E g 
requires a non-perturbative calculation of the correlation function of the boson field in the 
sine-Gordon field theory at finite temperature. Despite remarkably many exact results on the 
theory based on the integrability, calculation of the finite temperature correlation remains 
an unsolved problem. The ESR lineshape in Cu benzoate provides a set of rather precise 
experimental data for the finite temperature correlation function in the sine-Gordon field 
theory. It is hoped that future theoretical progress will enable us to compare theoretical 
non-perturbative result with the ESR data in the crossover temperature regime T ~ E g . 

In Ref. 11, the sharp resonance at very low temperature is considered to be the "antifer- 
romagnetic resonance," which reflects the Neel ordering due to the interchain interaction. 
In particular, they identified the appearance of the sharp resonance at very low temperature 
as the Neel transition. However, a recent /zSR experiment on Cu benzoate reveals 51 that a 



exp(-E g /T). 




(8.27) 




(8.28) 
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Neel ordering does not occur even down to 20mK. We believe that the evolution of the ESR 
lineshape is primarily explained within our one-dimensional theory taking the effective stag- 
gered field into account. On the other hand, we also note that the interpretation in Ref. 11 is 
not totally different from ours; the system has a long-range magnetic order in both theories. 
The difference is that the order is induced spontaneously due to the interchain interaction in 
Ref. 11 while it is forced externally by the staggered field in our picture. Spin-wave theory 
can also be applied to the externally ordered state; the resonance at very low temperature 
would be then identified with the "upper mode" E + (see eq. (3.12) of Ref. 7) which has a 
qualitatively similar dependence on h and H to eq. (8.27). (We thank H. Shiba for pointing 
this out.) Quantitatively, however, the sine-Gordon field theory is expected to work better 
for a small staggered field h. 



IX. ESR AT HIGHER TEMPERATURES 



In this paper, we have developed a field-theory approach to ESR in quantum spin chains. 
The field theory is a low-energy effective theory, and is only valid at low temperatures 
compared to the exchange coupling. Here we would like to consider ESR in the other 
extreme, namely the high temperature limit using Kubo-Tomita theory. We will also discuss 
the crossover between the low-temperature and high-temperature regime. 



A. Exchange anisotropy 




5 2 /y 



FIG. 6. Simplest scenario of the temperature dependence of the linewidth in the presence of 
an exchange anisotropy. The T-linear behavior at low temperature predicted by the field theory 
method crossovers smoothly to the constant ~ 5 2 / J predicted by the Kubo-Tomita theory at high 
temperature limit. The crossover takes place at T ~ J, which is the limit of the validity of the 
field theory approach. 



43 



For the exchange anisotropy (3.34) in a generic direction, the KT formula (2.12) has 
been applied to the linewidth in the literature. The result is 

5 2 

??oc 7 , (9.1) 

where we have ignored the direction dependence. It is difficult to discuss the intermediate 
temperature regime either with the existing theories or with our field theory approach. 
However, our result can be naturally related to the high-temperature limit, assuming a 
smooth crossover at temperature T ~ J, namely if the T-linear behavior (4. 8), (6. 27) is cut 
off at T ~ J as shown in Fig. 6. In fact, this simple scenario seems to agree with the 
experimental results 35-37,34,39,40 on CPC, KCUF3, CuGeC>3 and NaX^Os, which we think the 
exchange anisotropy (including the dipolar interaction) is the primary mechanism of the 
broadening. See Fig. 4 for some of the examples. 

We note that, while the low-temperature asymptotic behavior described by the field 
theory is universal, the crossover to the high-temperature regime is expected to be non- 
universal. The linewidth as a function of the temperature would depend, for example, on the 
next-nearest-neighbor interaction introduced additionally to the standard Hamiltonian (3.1). 



B. Staggered field 




FIG. 7. Simplest scenario of the temperature dependence of the linewidth in the pres- 
ence of a staggered field, interpolating the low-temperature field theory result Jh 2 /T 2 and the 
high-temperature result h 2 /J. 



No literature on the effect of a staggered field in ESR linewidth is known to us. The 
application of the KT formula (2.12) to the staggered field perturbation (8.1) is nevertheless 
straightforward, giving the linewidth 

h 2 

V oc —, (9.2) 

at the high temperature limit. This is again consistent with the low-temperature field theory 
result eq. (8.9), assuming a smooth crossover at T ~ J, as shown in Fig. 7. 
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C. Dzyaloshinskii-Moriya interaction 

While the effects of a DM interaction (3.32) has been discussed in the literature 33 ' 34,38 , we 
believe there is a rather serious problem with these previous treatments. A direct application 
of the KT formula (2.12), as was made previously, yields the width 

D 2 

where we have again ignored the angle dependence. 

On the other hand, a staggered DM interaction can be reduced to an exchange interaction 
5 ~ D 2 / J and a transverse staggered field h ~ DH/ J via an exact transformation discussed 
in Sec. Ill B 1. If we apply the KT formula after the transformation, we obtain 

D 2 H 2 

where we ignored constants of 0(1). This actually differs substantially from the result of 
the direct application (9.3). In a typical situation, D/J ~ 0.1 and H <C D < J so that 
i] ~ 0.01J from eq. (9.3) while rj ~ 1CT 4 J from eq. (9.4), which means a factor of 100 
difference. (Actually we have to know the numerical coefficient, which has been ignored so 
far, in order to discuss the absolute value of the width.) In addition, it is argued 28 that 
there exists an exchange anisotropy (before the transformation) which accompanies the DM 
interaction, and cancels the anisotropy coming from the DM interaction. The discrepancy 
would be even greater when this happens. 

Obviously, both results cannot be true at the same time (while they could be both 
wrong.) What we believe is that the latter approach eliminating the DM interaction first is 
appropriate, and the direct application of the KT formula to the DM interaction is incorrect. 
A possible reason why the direct application (9.3) fails is as follows. In the latter approach 
based on the transformation (9.4), the physical total spin operator S x,y is actually given 
by a sum of the total spin operator and the staggered spin operator of the model after the 
transformation. Thus, the physical absorption spectrum of ESR is also given by the sum of 
contributions from the uniform and staggered part: 

X ; hys (q = 0,u) ~ X "(q = 0,u)+ (y) X "(q = (9.5) 

where x" is the imaginary part of the dynamical susceptibility for the transformed model. 
The staggered part x"(q = 11 i u ) is already broad even in the absence of the anisotropic 
perturbation, and is further suppressed by the factor (D/J) 2 . Thus it would be practically 
indistinguishable from the background, especially in the high temperature regime. The 
main absorption due to the x"{q = 0)^) term is presumably Lorentzian with the width 
given by (9.4). According to this picture, the lineshape is not a single Lorentzian, although 
apparently it is. The direct application of the KT formula misses such a structure, and 
treats all the effects as if the lineshape is a single Lorentzian. This presumably leads to the 
incorrect result (9.3). 

An indirect evidence of our claim is that the elimination seems to work well in the field 
theory of ESR at low temperature. Assuming a smooth crossover at T ~ J, the latter 
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(9.3) 



(9.4) 



result (9.4) seems more plausible. In addition, a recent experiment 9 on a very good one- 
dimensional S — 1/2 Heisenberg antiferromagnet Pyrimidine Cu dinitrate strongly suggests 
that there is a staggered DM interaction along the chain, resulting in the field-induced 
gap similar to that observed in Cu benzoate. An analysis of various experimental data 
suggests 9 that the staggered DM interaction D ~ 0.14 J, where the exchange coupling in 
this compound is J ~ 36K. According to the direct approach (9.3), the linewidth at high 
temperature should be of the order of D 2 / J ~ 5000Oe. This might be too large to un- 
derstand the observed small linewidth ~ 20Oe at room temperature 9 , which is quite high 
compared to the exchange interaction J. On the other hand, if we use eq. (9.4) and H <C D, 
the estimate of the linewidth becomes to be of order of D 4 / J 3 ~ lOOOe, which is not too far 
from the experimental result. We note that we do not know the numerical coefficients and 
thus a conclusive quantitative discussion is difficult. In addition, the exchange anisotropy 
(before the elimination of the DM interaction), which is ignored in the above estimate, is not 
known precisely. Nevertheless, considering the significant difference, the observed linewidth 
in pyrimidine Cu dinitrate could serve as an experimental support for our claim that the 
direct treatment of the DM interaction is inappropriate. 

On the other hand, we do not understand at present how to deal with a uniform DM 
interaction along the chain. While it can be eliminated by a similar transformation as well, 
the result contains the magnetic field rotating in its direction along the chain. This is a 
rather unfamiliar problem which we do not know how to handle at present. 

D. High-temperature expansion of the linewidth 

In a series of papers, Yamada and collaborators studied the temperature dependence 
of ESR linewidth in one dimensional magnetic systems experimentally and theoretically. 
In the theoretical study, they discussed the temperature dependence by a high-temperature 
expansion of the KT formula. More precisely, they attempted a high-temperature expansion 
of (A4" 1 ") in the numerator of the KT formula, Eq. (2.12). 

They concluded that for an exchange anisotropy the linewidth increases as the tem- 
perature is lowered, while the tendency is the opposite for a (uniform or staggered) DM 
interaction. Based on this observation, they argued that the DM interaction should be dom- 
inant in several one-dimensional antiferromagnets which showed a decreasing linewidth at 
lower temperature. In some cases the DM interaction is forbidden according to the previ- 
ously identified crystal symmetry; they went on to the conclusion that the actual symmetry 
is lower than what had been believed, allowing the DM interaction. 

However, their argument is to be criticized on several grounds. First, the high- 
temperature expansion can not be trusted except for very high temperature. At T <C J 
our field theory approach should be more reliable, and it gives a rather opposite result to 
their claim. Second, they expand only the numerator (AA^) in the KT formula to the 
first order in 1/T, ignoring other possible contributions of order 1/T. It is not clear to us 
whether their scheme makes sense as a 1/T expansion of the linewidth. Third, perhaps most 
importantly, even in their framework of the calculation, the conclusion should be reversed 
because they apparently made a crucial sign mistake as we will show below. Finally, they 
apply the KT formula directly to the DM interaction; this is problematic as we have pointed 
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out. In any case, the sign problem persists whether the direct approach or the elimination 
approach is taken in dealing with the DM term. 

In the following, let us show that the sign should be reversed within the framework of 
Refs. 34, 38-40. We consider S — 1/2 Heisenberg antiferromagnetic chains with a small 
perturbation. First let us discuss the case of an anisotropy parallel to the applied field. The 
calculation for general anisotropy angle should be similar. This gives the commutator 

A=[H',S+] = 5Y,(S]S* +1 + S]Sf +1 ). (9.6) 
j 

The "numerator" of the KT formula is then given by 

(A4t) = 5 2 £<(^+i + S ]S+ +1 )(S k S z k+1 + S z k S k+1 )). (9.7) 

Considering the high-temperature limit we can ignore all but nearest-neighbor correlations. 
Thus we only consider the j = k terms in the double sum; 

W, S+][S-,H>]) ~ 5'Nl 1 - + ^ + (S*S* +1 )}, (9.8) 

where we have used identities for S — 1/2, such as SfSj — S? + 1/2 and N is the number 
of sites. In the limit of infinite temperature, the width is given by the first term which is 
5 2 / J as was already discussed. 

As the temperature is lowered from infinity, the leading correction is given by the second 
and third terms. The second term proportional to the magnetization (S z ) is negligible 
compared to the third term in our case H <C J. The third term represents the nearest- 
neighbor correlation effect, and should be proportional to —J/T at high temperature T. 
Note that we are dealing with an antiferromagnet, so that the nearest-neighbor correlation 
should be negative. Ignoring other possible sources of temperature dependence following 
Refs. 34, 38-40, the linewidth in the present case is given by 

V=j[a-bl + 0(£)], (9.9) 

with positive coefficients a, b for an antiferromagnet. Namely, the linewidth decreases at 
lower temperature contrary to the claims made in Refs. 34, 38-40; this is rather natural 
from the field theory results at low temperatures as discussed in Section IX A. It appears 
that, they took the nearest-neighbor correlation as positive, which is valid for a ferromagnet 52 
but not for an antiferromagnet. 

Now let us consider the transverse staggered field perturbation in the same framework. 
In this case, A = [H', S + ] = —h^2j(—iySj, which gives the "numerator" 

(AA*) = h 2 J2(-iy +k (S*S z k ). (9.10) 

In the high-temperature limit, we may ignore all the correlation functions other than the 
nearest neighbor one. This leads to the formula 
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(A4t) = h 2 N[\ - 2(S]S] +1 )] (9.11) 

Considering that the nearest-neighbor correlation is negative for an antiferromagnet, the 
linewidth is supposed to be given as 

V = j[a' + b'± + 0(f 2 )}, (9.12) 

where a' and b' are positive constants. Namely, the linewidth increases at lower temperature; 
again in a qualitative agreement with the field theory. 

For a staggered DM interaction, as we have discussed before, presumably we should 
first eliminate the DM interaction to reduce the problem to the exchange anisotropy and 
the transverse staggered field. In a typical situation H <C D <C J, the staggered field 
h ~ DH/ J and the anisotropy 5 ~ D 2 / J induced by the transformation satisfy h <C 5 <C J '. 
In this case, the linewidth would initially decrease by lowering the temperature, then start 
increasing below the crossover temperature where the staggered field becomes dominant. 
This was actually observed in Cu benzoate, as discussed already in Section VIII D. 

Finally, we consider a direct application of the KT formula to the DM interaction. Al- 
though we believe this is not an adequate approach, the claims in Refs. 34, 38-40 still suffers 
from the same sign problem even if we accept the direct approach. Now we have 

A= [H',S + ] =J2Dji(S*S+ +1 -S+S* +1 ), (9.13) 

3 

giving 

(AA ] ) ~ D 2 N[1 - (S*S* +1 )}, (9.14) 

ignoring other than next-nearest-neighbor correlation in the high-temperature limit. Because 
the nearest neighbor correlation function is negative in an antiferromagnet, we obtain 

V=^[a" + b"^ + 0(f 2 )], (9.15) 

where a" and b" are positive constants, implying the increasing linewidth at lower temper- 
atures. The error in Refs. 34, 38-40 is again apparently due to the identification of the 
nearest-neighbor correlation as positive. 



X. CONCLUSIONS 



In this paper, we have developed a new approach based on field theory to ESR in quantum 
spin chains. It is expected to be exact in the low-energy (low-temperature) limit, precisely 
where the traditional calculational methods on ESR become invalid. The weakly broken 
577(2) symmetry under an applied field, in the absence of an anisotropic perturbation, is 
represented by the SU (2) symmetric field theory and an anisotropic mapping between the 
physical spin operators and the corresponding field theory operators. 
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The formulation of the ESR in terms of Feynman-Dyson self-energy gives, at least in 
some simple cases, a microscopic derivation of the Lorentzian lineshape up to a possible 
smooth weak background. The spin diffusion picture 4 , which predicts a non- Lorentzian 
lineshape in one dimension, does not apply to the S — 1/2 antiferromagnetic chain at low 
temperature. The spin diffusion hypothesis does not hold in the present the spin 

correlation function is given explicitly using eq. (3.19). 

The width and shift are calculated perturbatively for a transverse staggered field pertur- 
bation and an exchange anisotropy parallel or perpendicular to the applied uniform field. 
They seem to explain many existing experimental data. Furthermore, the self-energy formu- 
lation can be used beyond the perturbation theory. In fact, in the presence of a staggered 
field, the perturbation theory breaks down at a low enough temperature. The ESR spectrum 
in the zero temperature limit is discussed with a non-perturbative treatment of the sine- 
Gordon field theory. This again seem to explain the experimentally observed ESR lineshape 
in Cu benzoate at very low temperature. 

While our field theory approach works only at low temperatures, we have also discussed 
a few aspects of ESR at higher temperatures. In particular, we have pointed out that a 
naive application of the standard Kubo-Tomita theory fails even in the high temperature 
limit, in the presence of a Dzyaloshinskii-Moriya interaction. 

We hope that the reader is convinced that ESR in a strongly interacting quantum sys- 
tem is quite an interesting problem from the theoretical point of view. It is also a useful 
experimental probe because a very precise spectrum can be obtained. 

Obviously, there remain many problems to be investigated in the future. Even in the 
simple quantum antiferromagnetic chain, the formulation of ESR in terms of self-energy of 
the boson field does not hold for generic types of anisotropic perturbations, because of the 
mixing of several operators. Extension of the self-energy formulation to the generic cases 
is an important open problem; presumably we have to consider perturbative expansion of 
correlation functions of the vertex operators (exponentials of the boson field) in a systematic 
way. Moreover, degrees of freedom other than spins (e.g. charge fluctuation, lattice vibra- 
tion etc.) will be relevant in some real materials. While the ESR in a three-dimensional 
magnet appears to be understood with the existing theory 2 ' 3 , we think that the problem 
should be reinvestigated with the modern understanding of many-body physics and critical 
phenomena. Naturally, the two-dimensional problem, which is expected to be more sensitive 
to the fluctuation effects, would also deserve consideration. We hope the present work will 
stimulate further theoretical and experimental studies on this fascinating subject. 

Note Added. After submitting the present paper, a paper by Choukroun, Richard and 
Stepanov was published 53 . They made a similar proposal to ours (Sec. IXC in the present 
paper) on the treatment of the Dzyaloshinskii-Moriya interaction at high temperature limit. 

ACKNOWLEDGMENTS 

It is a pleasure to thank Y. Ajiro, J. -P. Boucher, F. H. L. Essler, R. Feyerherm, Y. 
Maeda, S. Miyashita, Y. Natsume, H. Nojiri and H. Shiba for stimulating discussions and 
useful correspondences. In particular, we are grateful to D. M. Edwards for his suggestion on 
the content of Appendix and for allowing us to present it in this paper. We also thank A. A. 



49 



Zvyagin for pointing out Ref. 41 . This work has been supported in part by a Grant-in- Aid 
from MEXT of Japan and from NSERC of Canada. 



APPENDIX A: ALTERNATIVE DERIVATION OF THE MORI-KAWASAKI 

FORMULA 

In this appendix, we describe a simple alternative derivation of the MK for- 
mula (2. 7), (2. 10) suggested to us by David M. Edwards. It depends only on the assumption 
that the lineshape takes a single Lorentzian form, and appears much simpler than that in 
the original paper 2 . On the other hand, it does not answer the question why (and when) 
the lineshape takes the Lorentzian form. 

We consider ESR in a general spin system given by the Hamiltonian (2.6). Here and in 
the following, a spin operator without a site index is regarded as the total spin operator 
S a = J2j Sf- The equations of motion for S ± are given by 

= -iHS + + iA (Al) 

at 

^- = +iHS- - iA\ (A2) 
at 

where A= [H',S + ]. 

The ESR spectrum can be obtained from the Green's function of S ± . Let us relate this to 
the Green's function of A and A\ using the equations of motion. Using a partial integration 
and the equations of motion, 

PCX) 

g& s -(u,) = -i / e^([S+(t),S-(0)])dt 
Jo 

= -<[S+(0),S-(0)]> + -/ e**([—(t),S-(0)])dt 
UJ uj Jo dt 

= 2_{S^ + H g n +s {uj) _ Ie* 5 _ H . (A3) 

UJ UJ UJ 

Thus 

4-H = ^. (A4) 

UJ — H 

(Precisely speaking we should introduce the convergence factor so that uj — H is replaced 
by uj — H + ie with a positive infinitesimal e. Although we omit this for brevity, it can be 
recovered when necessary.) Performing similar steps, 



roc 

■ R (UJ 



AS-K") = -i e^{[A(t),S-(0)])dt 

J 

1 1 r°° d1~ 

-([A(0), S-(0)]> - - / e™([A(t), ^-(0)])dt 
uj uj Jo dt 



(L4(0),£-(0)]) H 1 

+ —yAS- \V) - -^AAt ( w ) > ( A5 ) 



+ 

uj uj ~" 1 ' UJ 
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where we used the relation 



j t ([A(t),S-(0)}) = -([A(t)^(0)}) (A6) 

which holds because the Green's function depends only on the difference of two time argu- 
ments. Thus 

&H J» S 'W'-^'"' (A7) 

uo — H 



Combining eqs. (A4) and (A7), we obtain 

(A8) 



c r ( v 2(S*) -([A(0),S-(0)]) + Q%(u) 



This should be an exact relation between the full Green's functions (in which the effect of 
the perturbation W is fully taken into account.) When 7i' = 0, we recover the simple result 
gfj+ s -(u) = 2(S')/(u-H). 

Now let us assume the perturbation H! is small, and the ESR lineshape is Lorentzian. 
Namely, we assume that Qg +s -{uj) is given by eq. (2.11) where £ is a smooth function of uj. 
Near the resonance u ~ if, E may be regarded as a constant. ReS and — ImS gives the shift 
and width of the resonance, respectively. We assume that S can be expanded perturbatively 
in W. 

Comparing eqs. (2.11) and (A8), we obtain, in the lowest order of perturbation theory, 
v _ -([A(0),S-(0)])+S 5a^ = H) ^ (A9) 



2(S Z ) 

Here we note that ([^4(0), S'~(0)]) is purely real since [A, S~\ is Hermitean. This gives 

V = ^ Im ^t(^ = #), ( A1 °) 



Ao; = 



1 



2(S Z ) 



-([A,S-]) + Reg^(u = H)}. (All) 



For a small field H, the denominator 2(S' 2 ) can be written as 2x u H where \ u is the uniform 
susceptibility. We also note that the first term in the shift — ([A, S~]) / (2(S Z )) was derived 
previously by Kanamori and Tachiki 47 , and by Nagata and Tazuke 48 . However, their theory 
did not incorporate the dynamical effects represented by Gj^- 

So far, we have defined the expectation value and the Green's functions with respect to 
the full Hamiltonian 7i = Tio + Tiz + H! . However, the present result is only valid in the 
leading order. Since the Green's function of A above already contains the factor A 2 (A is 
the small parameter that characterizes the perturbation H'), we may replace G%^ with the 
unperturbed Green's function G^t- This gives the formulae in eqs. (2.7) and (2.10). We 
note that, in general, ([A, S~]) must still be evaluated in the presence of H' because [A, S~] 
is only first order in A. 
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